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ABSTRACT 


This  study  was  initiated  to  determine  the  effects  of  antenna  element 
weighting  (amplitude  control)  and  spacing  on  the  signal-to-noise 
ratio  of  wideband  systems.  The  initial  analysis  is  directed  toward 
the  linear  array  operating  in  an  environment  of  white  noise,  uniformly 
distributed  in  space.  These  confinements  are  later  removed  to 
illustrate  application  to  planar  arrays,  and  non-uniform  noise  con¬ 
ditions. 

It  is  shown  that  the  signal-to-noise  ratio  of  an  antenna  system  operating 
in  a  uniform  noise  environment  is  optimized  when  the  array  directivity 
is  maximum.  Optimization  involves  the  application  of  a  specific 
illumination  function  which  can  be  achieved  by  element  weighting,  ele¬ 
ment  spacing  ora  combination. 

The  application  of  wideband  signals  modifies  the  array  pattern  in  the 
time  space  domain.  Although  the  nature  of  the  modification  is  de¬ 
pendent  on  waveform,  the  sidelobe  displacement  in  time  and  space 
is  a  function  of  the  number  of  elements  comprising  the  array. 

An  iterative  technique  of  mathematical  analysis  is  developed  which 
permits  solutions  to  large  array  problems  where  matrix  inversion 
is  difficult.  Several  examples  of  possible  applications  of  the  analysis 
are  presented  including  a  spectral  detection  technique,  an  approach 
to  secure  communications,  and  a  concept  of  array  adaptation  to  a 
non-uniform  noise  environment. 


EVALUATION 


Much  study  has  been  devoted  to  optimization  of  element  spacing 
and  element  weighting  for  narrow  band  arrays,  but  comparatively 
little  equivalent  study  has  been  devoted  to  the  wideband  array. 

The  study  herein  hai  looked  at  the  wideband  array  and  has 
further  developed  the  theory  of  wideband  arrays  by  considp; ing 
not  only  radiation  patterns  tut  also  the  processed  array  signal. 

Work  in  the  area  of  wideband  arrays  will  continue  to  put  the 
theory  and  techniques  of  the  wideband  array  on  the  same  firm 
basis  as  the  narrow  band  array. 


/JOSEPH  LOVECCHIO 


Project  Engineer 
EMATS/330-3685 


.  //  <1.  .  \ 

THOMAS  MAGGIO  // 

Chief,  Signal  Processing  Section 
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SECTION  I 


INTRODUCTION 


A  wideband  array  is  defined  as  one  in  which  the  excitation  consists  of  some¬ 
thing  other  than  a  pure  continuous  wave  sinusoid.  In  general,  the  excitation  may  take 
any  form,  from  a  pair  of  discrete  frequency  components  to  a  continuous  spectrum  de¬ 
fined  in  terms  of  its  waveform. 

The  purpose  of  this  investigation  is  to  determine  the  response  characteristics 
of  array  antennas  excited  by  wideband  signals.  The  array  and  the  applied  excitation 
are  considered  variables  in  an  attempt  to  optimize  the  signal-to -noise  ratio  of  the 
processed  array  output.  The  primary  array  variables  are  the  separation  between  ele¬ 
ments  (spacing) ,  and  the  relative  amplitude  coefficients  at  the  element  terminals 
(weighting) . 

In  order  to  provide  a  hasis  for  comparison,  a  set  of  initial  conditions  is  as¬ 
sumed  which  confines  the  extent,  but  not  the  scope,  of  the  investigation. 

a.  Linear  arrays  are  to  be  considered,  although  no  restraint  is  placed 
on  size,  number  of  elements,  or  spacing. 

b.  Only  fixed  beam  arrays  with  a  broadside  beam  are  considered. 

Analysis  must  be  carried  over  sufficient  coverage  to  detect  the 
presence  of  any  grating  lobes. 

c.  A  point  signal  source  is  assumed  to  be  aligned  with  the  beam  maximum. 

d.  The  source  of  noise  is  assumed  to  be  "white1 '  (evenly  distributed  in 
frequency)  said  uniformly  distributed  over  all  space. 

e.  Control  of  the  signal-to -noise  ratio  is  effected  by  manipulation  of  the 
element  spacing  and  element  weighting. 

f.  Wideband  signals  are  applied  to  determine  the  effects  of  bandwidth, 
waveform ,  and  spectral  content. 

The  major  portion  of  this  report  concerns  array  response  within  these  condi¬ 
tions;  however,  the  linear  array  restriction  (a) ,  and  the  uniform  noise  restriction,  (d) , 
are  removed  to  illustrate  the  general  applicability  of  the  results. 

Evaluation  of  the  array  modification  techniques  is  achieved  in  comparison  with 
a  standard  model  array.  The  standard  model  is  defined  as  a  linear  array  of  N  uniformly 
illuminated  isotropic  elements  located  on  half-wavelength  centers  at  a  fixed  frequency. 

The  gain,  beamwidth,  sidelobe  level,  and  signal -to -noise  ratio  (in  white  noise)  are  fixed 
by  N  for  this  standard  array. 

In  an  attempt  to  provide  a  coherent  presentation  of  the  results  obtained  during 
the  investigation,  a  step-by-step  reporting  procedure  is  incorporated.  This  procedure 
begins  with  an  analysis  of  narrow  band  linear  arrays,  with  control  provided  through 
amplitude  weighting.  Element  displacement,  to  effect  a  "space  taper"  is  also  discussed. 
The  next  step  is  to  consider  the  effect  ot  noise  on  array  performance  and  to  establish 
signal-to -noise  relationships  utilizing  the  array  parameters.  The  array  conditions  thus 
developed  are  subjected  to  the  effect  of  wideband  signals  so  that  a  composite  picture  of 
the  wideband  optimized  array  is  produced.  The  results  of  this  investigation  have  provided 
several  features  of  the  wideband  array  which  suggest  application  in  practical  array  sys¬ 
tems.  These  app'ications  are  discussed  in  a  separate  section. 
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SECTION  II 


SUMMARY 


The  principal  objective  of  the  Wideband  Element  Spacing  and  Weighting  inves¬ 
tigation  is  to  determine  the  procedure  whereby  the  signal-to-noise  ratio  of  a  wideband 
antenna  array  can  be  optimized  by  means  of  antenna  element  amplitude  control  and 
spacing.  A  subsequent  objective  is  to  determine  the  array  response  as  these  parameters, 
along  with  bandwidth,  are  varied  over  a  prescribed  range  of  values. 

Some  of  the  results  of  this  investigation  can  be  expressed  as  statements  which 
are  substantiated  in  the  body  of  the  report. 

a.  Maximum  signal-to-noise  ratio  is  synonomous  with  maximum  array 
directivity  when  external  noise  is  uniformly  distributed  in  frequency 
and  space. 

b.  Maximum  signal-to-noise  ratio  in  narrow  band  rays  can  be  specified 
by  an  optimum  illumination  function  for  a  spec  .ied  sidelobe  level. 

c.  The  optimum  illumination  can  be  applied  as  an  amplitude  weighting 
function,  a  space  taper  function,  or  a  combination  function. 

d.  An  improvement  in  signal-to-noise  ratio,  over  the  uniform  model, 
can  be  obtained  by  applying  the  super -gain  effect  if  the  resulting 
array  parameters  can  be  established. 

e.  The  application  of  bandwidth  modifies  the  array  function  in  the  time- 
space  domain. 

f .  Beamwidth  can  be  incorporated  with  amplitude  or  space  weighting  to 
produce  specifically  shaped  array  patterns. 

g.  Increased  bandwidth  reduces  the  optimum  array  signal-to-noise  ratio, 
but  permits  optimization  within  a  confined  aperture. 

h.  Optimized  amplitude  weighting  deviates  from  uniform  only  slightly, 
even  in  wideband  arrays. 

i.  An  iterative  technique  can  be  used  to  solve  large  array  problems  where 
matrix  inversion  becomes  cumbersome. 

j.  The  optimized  array  response  of  the  planar  array  is  similar  to  its 
counterpart  linear  array. 

k.  Waveform  distortion  is  prevalent  in  the  sidelobe  response  of  wide¬ 
band  arrays. 

l.  The  array  signal-to-noise  ratio  can  be  optimized  for  non-uniform 
noise  conditions  by  amplitude  control,  resulting  in  a  modified  array 
pattern. 

The  design  parameters,  which  control  the  performance  of  a  linear  array,  are 
discussed  in  Section  III.  Most  of  the  information  is  presented  in  graphical  form  so  that 
an  estimate  of  array  performance  can  be  obtained  without  computation.  An  example  is 
provided  to  illustrate  the  process. 
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The  effects  of  noise  on  array  performance  is  discussed  in  Section  IV.  A 
general  noise  analysis  is  followed  by  signal -to -noise  ratio  considerations  and  several 
optimization  techniques  based  on  element  amplitude  weighting  and  spacing. 

Bandwidth  is  considered  in  Section  V.  The  effect  of  bandwidth  on  the  radia¬ 
tion  patterns  in  the  time-space  domain  is  illustrated  by  a  three-dimensional  model 
and  a  series  of  contour  diagrams.  The  effect  of  bandwidth  on  array  signal-to-noise 
ratio  is  presented  in  graphical  form.  An  application  of  the  space -frequency  equivalence 
theorem  applied  to  a  space  tapered  array  is  presented  as  an  example  of  utilizing  band¬ 
width  in  array  design. 

Several  possible  applications  of  the  principles  evolving  from  the  analytical 
investigation  are  discussed  in  Section  VI.  These  applications  are  intended  to  illustrate 
the  manner  in  which  the  information  developed  during  this  project  can  be  applied  in 
future  development  effort. 

Several  computational  techniques  have  been  developed  during  the  course  of 
investigation.  Since  these  methods  are  general  in  nature,  they  are  discussed  sepa¬ 
rately  in  Section  VII. 

Considerable  information  and  several  mathematical  formulations  have  evolved 
which  supplement  the  data  and  statements  in  the  body  of  this  report.  These  are  col¬ 
lected  as  appendices  to  the  report. 

A  list  of  references  and  a  bibliography  are  provided  to  permit  further  evalu¬ 
ation  of  the  general  topics  covered  in  this  report. 
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SECTION 


III 


DESIGN  PARAMETERS  IN  THE  LINEAR 
ARRAY  --  NARROW  BAND  ARRAYS 


An  application  of  wideband  element  spacing  and  weighting  techniques  to  a 
specific  antenna  array  problem,  requires  information  concerning  the  performance 
of  the  array  for  the  specified  parameters.  In  this  section,  the  parameters  which 
determine  the  performance  of  an  array  are  examined  with  emphasis  on  their  inter¬ 
dependence.  The  narrow  band  restriction  is  applied  to  permit  an  independent  ex¬ 
ploration  of  array  parameters,  so  that  the  application  of  bandwidth  (Section  V)  can 
be  considered  in  proper  perspective. 

An  antenna  array  is  usually  specified  in  terms  of  its  gain,  beam  width,  and 
sidelobe  level.  There  are  other  parameters  which  may  be  specified  such  as  polar¬ 
ization,  available  aperture,  method  of  excitation,  input  impedance,  power  handling 
capability  and  form  factor,  but  these  are  either  related  to  the  three  primary  para¬ 
meters  or  can  be  considered  independently.  One  parameter  which  is  an  overriding 
function  in  any  antenna  design  is  efficiency.  Efficiency  includes  all  factors  which 
represent  a  departure  from  maximized  directivity,  which  ir.  this  report  is  defined 
as  a  lossless  uniformly  illuminated  array. 

In  a  practical  design,  the  array  efficiency  affects  all  of  the  above  para¬ 
meters  eithe-  directly  or  indirectly.  For  purposes  of  discussion  in  this  section 
only  the  illumination  efficiency  will  be  considered  in  analysis  sii.ce  system  losses 
can  be  applied  on  a  discrete  basis. 

1.  ARRAY  DESIGN 

The  gain  of  an  array  is  always  a  function  of  the  aperture  available.  Within 
the  specified  aperture  .beam  width  and  sidelobe  level  can  be  controlled  by  means  of 
the  applied  illumination  function. 

The  radiation  pattern  of  an  array  (which  provides  beam  width  and  sidelobe 
level)  can  be  synthesized  from  a  known  or  assumed  illumination  function  and  the 
geometrical  properties  of  the  array.  The  illumination  function  can  be  described 
as  a  continuous  function  or  as  a  pattern  of  discrete  weighting  factors.  The  illum¬ 
ination  function  can  be  implemented  as  an  amplitude  weighting  on  individual  elements 
(amplitude  taper),  as  a  non-uniform  distribution  of  equal  amplitude  elements  w'ith- 
in  the  array  (space  or  density  taper)  or  as  a  combination  of  these  techniques. 

The  details  of  antenna  array  synthesis  have  been  adequately  covered  in 
the  available  literature*  and  will  not  be  repeated  hei’e.  An  analytical  technique 
which  is  applicable  to  array  analysis  and  has  not  been  previously  reported  appears 
in  Section  VII.  The  purpose  of  this  discussion  is  to  illustrate  the  application  of 
specific  illumination  functions  to  the  general  form  of  the  radiation  pattern.  The 
selected  functions  are  briefly  described  below  and  nre  tabulated  in  numerical  form 
in  Tables  1  and  2  for  continuous  and  discrete  weighting  respectively. 


*An  excellent  survey  of  the  theory  of  antenna  arrays  is  contained  in  Reference  ( 1). 
NOTE:  Reference  numbers  noted  refer  to  references  listed  in  Section  VIII. 
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TABLE  1.  CONTINUOUS  APERTITIE  WEIGHTING 
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TABLE  1.  CONTINUOUS  APERTURE  WEIGHTING 
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1.  1  UNIFORM  WEIGHTING 


This  condition  represents  maximum  directiviiy  and  is  assigned  an  illumin¬ 
ation  efficiency  of  100  percent.  Uniform  illumination  can  he  effected  on  a  continuous 
or  discrete  basis.  In  either  case  the  first  sidelobe  is  approximately  13  DB  below 
the  major  lobe.  The  difference  in  the  two  forms  of  uniform  illumination  is  observed 
in  the  far  out  sidelobe  pattern  where  the  continuous  illumination  exhibits  a  lower 
level.  This  difference  is  accentuated  as  the  total  number  of  elements  is  reduced. 

A  discrete  array  of  only  four  elements,  foi  example,  will  produce  a  sidelobe  only 
11  DB  below  the  beam  maximum  rather  than  the  characteristic  13  DB  of  larger 
arrays.  Kraus,  Reference  (  2 )  has  recorded  the  space  patterns  for  discrete  uniform 
illumination  in  an  array  of  N  elements  where  N  is  any  integer  from  1  through  24. 

1.  2  HAMMING  FUNCTION 

Named  for  R.W.  Hamming  -  Bell  Laboratories,  Reference  (3  ),  the  Ham¬ 
ming  Function  like  uniform  weighting  can  be  applied  on  a  uniform  or  discrete  basis. 
The  function  is  essentially  a  cosine  variation  in  a  form  which  represents  a  compro¬ 
mise  between  major  lobe  broadening,  sidelobe  suppression,  illumination  efficiency 
and  ease  of  implementation. 

1.  3  COSIN En  ILLUMINATION 

The  cosine  function  is  a  basic  form  of  illumination.  Sidelobe  control  is 
exercised  by  variation  of  the  exponent  "n".  Characteristics  of  cosine11  illumination 
are  significant  beam  broadening  and  low  efficiency.  These  undesirable  character¬ 
istics  can  be  improved  by  placing  the  function  on  a  pedestal,  i.  e.  restricting  the 
edge  term  to  a  non-zero  value.  The  cosine11  on  a  pedestal  resembles  the  Hamming 
function. 

1.4  MODIFIED  SIN  X/X 

The  name  applied  to  the  resulting  radiation  pattern  rather  than  the  illum¬ 
ination  function.  This  method  of  array  excitation  was  suggested  by  Taylor, 

Reference  (  4  ),  as  a  means  of  achieving  the  rapid  sidelobe  decay  of  the  uniform 
array  without  the  high  first  sidelobe  level. 

1.  5  TRIANGULAR  FUNCTION 

This  is  included  primarily  for  comparison  purposes.  The  triangular 
illumine  ion  function  is  seldom  used  in  practice  although  the  sidelobe  level  can  be 
controlled  by  establishing  the  magnitude  of  the  edge  illumination.  For  low  side- 
lobes,  a  broad  beam  with  low  efficiency  results. 

1.  6  TCHEBZCHEFF  DISTRIBUTION 

A  technique  developed  by  Doiph,  Reference  (5  ),  of  applying  the  Tchebzcheff 
polynominals  to  antenna  design  results  in  uniform  sidelobes  with  a  minimum  beam 
width.  This  form  of  excitation  has  found  extensive  application  despite  its  low  illu¬ 
mination  efficiency. 

1.7  TAYLOR'S  DISTRIBUTION,  (Reference  (  6 )) 

An  adaptation  of  the  Dolph-Tchebzcheff  form  which  controls 
a  specified  numbr^  of  sidelobes  to  a  uniform  level  while  the  remainder  decay  with 
increased  displacer  nt  uom  the  major  lobe  is  a  distribution  form  which  has  found 
extensive  use  in  radar  antennas. 
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The  above  described  functions  serve  to  illustrate  the  pattern  control  which 
is  available  through  amay  illumination.  There  are  many  other  illumination  func¬ 
tions  which  can  be  applied  including  modifications  of  the  above.  The  selection 
must  be  based  on  the  desired  radiation  characteristics  which  are  governed  by 
application. 

2.  GRAPHICAL  ANALYSIS 

The  analytical  forms  presented  in  Tables  1  and  2  are  complete  and 
permit  direct  numerical  analysis  for  the  functions  tabulated.  For  a  complete  array 
design,  the  specified  parameters  should  be  applied  through  the  selected  illumination 
function  io  produce  the  design  array  function.  For  purposes  of  preliminary  analy¬ 
sis,  however,  it  is  often  desirable  to  determine  the  approximate  characteristics 
before  proceeding  with  a  detailed  design.  This  is  particularly  true  when  other 
factors  are  to  be  considered,  such  as  those  presented  in  subsequent  sections  of 
the  report. 

The  pattern  characteristics  resulting  from  applied  illumination  functions 
can  be  displayed  in  graphic  and  nomographic  form  provided  a  reference  is  estab¬ 
lished.  The  reference  used  in  this  section  is  the  uniformly  illuminated  array. 

The  calculated  array  function  is  illustrated  in  Figure  1  along  with  the  bidelobe 
envelopes  for  other  illumination  functions. 

The  illumination  efficiency  of  the  uniformly  illuminated  array  is  taken  to 
be  100  percent.  The  sidelobe  behavior  is  shown  in  Figure  1  and  the  beam  width 
as  a  function  of  array  length  as  can  be  determined  from  Figure  2  .  The  rela¬ 

tive  illumination  efficiency  for  other  illumination  functions  can  be  determined  from 
Figure  3  as  a  function  of  design  sidelobe  level.  The  half-power  beam  width  of 
the  array  can  be  determined  from  Figure  4  which  gives  the  beam  width  multi¬ 
plication  factor  as  a  function  of  sidelobe  level  for  several  illumination  functions. 

The  directive  gain  of  the  linear  array  is  determined  by  the  number  of 
elements  comprising  the  array,  the  illumination  efficieney,  and  the  gain  of  the 
individual  elements.  This  relationship  is  presented  as  a  nomograph  in  Figure  5 

Figures  x  through  5  can  be  used  to  provide  a  rapid  approximate  solu¬ 
tion  to  various  array  design  problems,  and  provide  a  visual  display  of  parameters. 
To  illustrate  the  usefulness  of  the  graphical  analysis,  consider  a  sample  problem. 

a.  Illustrative  Example 

It  is  desired  to  determine  the  gain  and  beam  width  of  a  70  element 
linear  array  with  half  wavelength  element  spacing  designed  to  provide 
a  maximum  sidelobe  level  of  25  DB  below  the  major  lobe.  Assume 
that  sidelobe  control  can  be  achieved  through  application  of  the  Ham¬ 
ming  illumination  function.  The  array  element  gain  is  8  DB. 

b.  Procedure 

(1)  Determine  the  array  length  by  multiplying  the  quantity  of  elements 
(70)  by  the  element  spacing  (A/2).  L  =  35A 

(2)  From  Figure  2  determine  the  half-power  beam  width  for  a 

35A  linear  uniformly  illuminated  array.  Beam  width  =1.4  degrees 
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Illumination  Functions 


HALF  POWER  BEAMWIDTH  (DEGREES) 


TOTAL  ARRAY  LENGTH  (WAVELENGTHS) 


Figure  2.  Beamwidth  of  Uniform  Linear  Array  (Beam  Normal  to  Array  Axis) 
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Figure  3.  Illumination  Efficiency  in  Amplitude  Tapered  Linear  Array 
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BEAMWIDTH  MULTIPLICATION  FACTOR 


Figure  4.  Beamwidth  Relative  to  Uniformly  Illuminated  Linear  Array 
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NUMBER  OF  PIVOT  ARRAY  ILLUMINATION  ELEMENT 

ELEMENTS  LINE  GAIN  EFFICIENCY  GAIN 

N  G(DB)  r,  ge(DB) 


PROCEDURE: 

1.  CONNECT  N  AND  i\  WITH  A  STRAIGHT  LINE  WHICH  INTERSECTS  THE 
PIVOT  LINE. 

2.  CONNECT  THE  PIVOT  LINE  INTERSECTION  WITH  THE  ELEMENT  GAIN  AND 
READ  ARRAY  GAIN  FROM  THE  INTERSECTION  WITH  THE  G  SCALE. 


Figure  5.  Linear  Array  Nornograph 


(3)  From  Figure  4  determine  the  beam  width  multiplication  factor 
for  a  Hamming  function  of  25  DB  sidelobe  level  design.  Multi¬ 
plier  =  1.  19 

(4)  The  product  of  steps  (2)  and  (3)  gives  the  approximate  half-power 
beam  width.  BW  (3  DB)  =  1.  68  degrees 

(5)  From  Figure  3  determine  the  relative  efficiency  for  the  25  DB 
sidelobe  Hamming  function,  r,  =  90% 

(6)  Using  the  nomograph  (Figure  5  ),  connect  75  on  the  N  scale 

with  90  percent  on  the  rj  scale. 

(7)  Mark  the  point  where  this  line  intersects  the  pivot  line  and  con¬ 
nect  this  point  to  8  on  the  gg  scale. 

(8)  Read  the  approximate  array  gain  from  the  G  scale.  G  =  26.  3  DB 

The  general  shape  of  the  radiation  pattern  in  the  sidelobe  reg’on  can  be 
seen  in  Figure  1  .  The  Hamming  function  sidelobe  envelope  is  for  a  sidelobe 
maximum  of  23  DB,  but  the  general  distribution  is  indicative  of  other  design  levels. 


The  parameter  ge  —  element  gain,  requires  a  word  of  explanation  in  appli¬ 
cation  to  linear  arrays-  In  a  two  dimensional  planar  array,  the  individual  radiating 
element  is  confined  to  a  physical  area  which  confines  the  effective  area;  hence  the 
gain.  In  the  linear  array,  confinement  is  limited  to  one  dimension  only,  therefore 
it  is  possible  to  utilize  high  gain  elements  in  the  linear  array.  An  approximation 
of  the  array  element  gain  can  be  established  from  a  knowledge  of  a  separation 
between  elements  and  the  half-power  beam  width  in  the  plane  normal  to  the  array. 
Figure  6  illustrates  this  relationship.  Using  the  above  example  where  the 
element  separation  was  0.  5A  and  the  element  gain  given  as  8  DB,  the  half  power 
beam  width  of  the  array  in  the  orthogonal  plane  can  be  found  from  Figure  6  to 
be  60  degrees.  The  depicted  relationship  is  only  an  approximation,  since  the  gain 
of  any  antenna  element  is  affected  by  factors  other  than  beam  width  and  element 
separation. 

3.  SPACE  TAPER  ARRAYS 

The  graphical  method  of  array  analysis  outlined  above  applies  only  to  the 
amplitude  tapered  array.  It  is  also  possible  to  effect  a  taper  function  by  the  dis¬ 
placement  of  elements  along  the  line  of  the  array  while  the  magnitude  of  excitation 
is  uniform.  In  general  this  procedure  consists  of  estimating  the  amplitude  or 
spacing  coefficients  necessary  to  produce  a  desired  pattern  and  calculating  correc¬ 
tion  factors  from  the  observed  pattern  necessary  to  produce  the  desired  result. 

The  calculated  patterns  shown  in  Figures  7  and  8  illustrate  the 
results  of  the  iterative  process  for  a  24  element  array  designed  for  24  DB  side- 
lobes.  The  two  figures  repsent  the  amplitude  and  space  tapered  cases  respective¬ 
ly.  Comparison  of  these  two  methods  of  pattern  control  reveals  the  following 
characteristics: 

a.  The  space  tapered  pattern  is  more  sensitive  to  small  deviations  in 
displacement,  than  the  amplitude  tapered  pattern  is  to  small  variations 
in  amplitude. 

b.  An  exact  sidelobe  level  is  more  easily  obtained  with  the  amplitude 
taper. 
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HALF  POWER  BEAMWIDTH  NORMAL  TO  ARRAY  (DEGREES) 
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iercd  Array  of  24  Elements  Sidelobe  Control  in  /u/<l.  00 


c.  The  beam  width  of  the  amplitude  tapered  array  increases  as  a 
function  of  s 'delobe  level  more  than  the  space  tapered  array  The 
beam  width  in  the  space  tapered  array  is  only  slightly  wider  than 
the  uniformly  spaced,  uniformly  illuminated  array. 

d.  The  illumination  efficiency  for  the  space  tapered  array  is  lower  than 
for  the  amplitude  tapered  array  with  a  common  sidelobe  level. 

e  The  element  gain  in  the  space  tapered  array  is  not  constant  since  the 
spacing  between  elements  is  not  constant.  Except  in  severely  tapered 
arrays  this  is  not  a  serious  factor  in  array  response. 

Although  the  graphical  method  of  array  evaluation  does  not  apply  to  the 
space  tapered  array,  a  fair  estimate  for  performance  can  be  obtained  by  assuming 
an  equivalent  amplitude  taper  condition  and  applying  the  above  difference  condi¬ 
tions. 
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SECTION  IV 


SIGNAL-NOISE  CONSIDERATIONS  IN  ARRAY  ANTENNAS 


The  procedure  for  antenna  array  analysis  described  in  the  preceding  section  con¬ 
sidered  only  the  signal  aspect  of  antenna  performance.  The  antenna  radiation  pattern 
is  an  expression  of  the  relative  signal  strength  (with  respect  to  the  beam  maximum)as  a 
function  of  observation  angle.  The  absolute  value  of  signal  available  for  use  in  a  detec¬ 
tion  circuit  is  dependent  on  a  number  of  other  factors  associated  with  the  antenna 
system.  These  include:  (1)  the  spatial  power  density  at  the  array  aperture,  (2)  the 
impedance  of  the  antenna  and  its  associated  components,  and  (3)  the  loss  or  gain  of 
any  subsequent  components.  Reference  (7). 

In  addition  to  the  signal  available  for  detection,  the  noise  assoc  atei  with  the 
system  and  its  environment  must  be  considered  since  the  practical  usefulness  of  any 
detector  is  its  ability  to  separate  signal  from  the  input  noise.  The  figure  cf  merit 
which  is  genially  used  in  analyzing  circuits  with  noise  is  the  signal -to-noise  ratio 
(S/N). 

Signal  and  noise  parameters  are  considered  in  this  section  as  they  relate  to 
narrow  band  array  antennas.  The  bandwidth  restriction  is  temporary  and  is  removed 
in  Section  V.  The  object  of  this  analysis  is  to  determine  what  control  of  the  overall 
signal-to-noise  ratio  can  be  obtained  through  manipulation  of  element  weighting 
(amplitude  control)  and  spacing  (density  control). 

1  SIGNAL  POWER  IN  THE  ANTENNA  ARRAY 

For  purposes  of  analysis  the  antenna  array  can  be  represented  as  a  "black  box" 
with  N  inputs  and  a  single  output.  The  concept  is  illustrated  in  Figure  9  .  The 
antenna  signal  response  can  be  defined  as  value  of  signal  power  (W)  which  can  be 
dissipated  in  a  terminating  impedance  (Zt)  for  an  applied  constant  fixed  field  density 

<$o> 


Since  the  antenna  response  is  a  function  of  angle,  it  is  appropriate  to  utilize  the 
concept  of  a  gain  function  which  is  defined  in  terms  of  the  array  operating  as  a  re¬ 
ceiving  antenna.  Assuming  reciprocity,  the  expressions  are  also  valid  for  the  trans¬ 
mitting  case.  The  definition  of  gain  is  the  ratio  of  the  power  available  for  dissipation 
in  a  toad  resistor  (Rt),  to  the  total  power  available  to  tip  antenna.  In  antenna  work 
it  is  convenient  to  express  this  ratio  as  a  relative  ratio  using  an  isotropic  source  as 
the  reference  (7). 


The 

where: 


gain  of  any  antenna  relative  to  an  isotropic  source  is  then  G  =  — 

W' 


W  =  power  delivered  to  load 

W1  =  power  delivered  to  load  if  antenna  is  an  isotropic  source. 


(1) 
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WAVE 

V  *  VOLTAGE  INDUCED  BY  THE  INCOMING  WAVE 
Za  =  ANTENNA  IMPEDANCE  »  R*  +  ptz 

I  =  CURRENT  FLOWING  AS  A  RESULT  OF  THE  INDUCED  VOLTAGE  V 
Zt  =  Rt  +  pet  (TERMINATING  IMPEDANCE) 

1  =  Za  +  Zt 

Figure  10.  Equivalent  Circuit 
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The  absolute  gain  of  any  antenna  is  a  function  of  the  effective  area  of  the  antenna 
aperture.  But  the  effective  area  of  an  isotropic  source  is  \2/4tt  therefore  the  power 
received  by  the  unknown  antenna  relative  to  that  which  would  have  been  received  by 
an  isotropic  source  is  related  to  the  gain  of  the  antenna  relative  to  the  gain  of  the  iso¬ 
tropic  source  (unity),  assuming  a  constant  fixed  field  intensity  . 

G  oc  4?r  W  (2) 

In  the  simple  aperture  case  the  antenna  can  be  represented  by  the  equivalent  circuit 
of  Figure  10. 

Considering  a  lossy  antenna  R  =  R  +  R  where  R  is  the  antenna  radiation  resist- 

2L  r  6  r 

ance  and  Rg  is  the  loss  in  the  antenna  and  circuit. 

The  total  pov/er  available  to  the  antenna  from  the  incoming  wave  is 


but  the  power  available  for  dissipation  in  the  load  impedance  (Z^)  is 


if  P  =  uW 


2  P\ 
W  =  1  Rt  =  "Jr 


Pi 

W 


v : 

R. 


l2<Za+Zt>2 

R, 


W-I  <za  +  Zt> 
2 

u  R. 


where  u  is  an  efficiency  factor.  (NOTE:  Kraus,  Reference  (7),  refers  to  u  as  a 
scattering  coefficient  in  the  lossless  antenna.) 


Substituting  for  Z  and  Z, 

cl  t 


w  = 


(R  +  R  +  R.)2  +  (X  +  XJ2 
'  r  e  v  v  a  r 


(3) 


2 

The  maximum  power  which  can  be  dissipated  in  the  load  is  I  Rt.  This  condition 
exists  only  when  X  =  X  ,  R  =0,  and  R  =  R  which,  from  equation  3,  renders 

3.  t  g  t  r 


u  =  2.  That  is  to  say  that  under  optimum  receiving  conditions  (conjugate  matching) 
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only  half  of  the  available  power  can  be  absorbed  in  the  load.  The  remainder  is 

scattered  from  the  antenna  (absorbed  in  R  ).  Under  these  conditions,  the  maximum 

r 

value  of  incident  power  which  can  be  absorbed  in  the  load  is 


<V 


R  + 


Rt> 


+  (X„  +  xY 


(4) 


Extending  this  simple  analysis  to  the  array  depicted  in  Figure  b  we  can  consider 
each  element  as  an  incremental  collector  as  shown  in  Figure  11. 


Following  the  preceding  analysis,  the  power  available  in  terminating  impedance 
th 


of  the  n  element  is 
i2 


W  =  7#-  (R  +  R  +  R,  )2  +  (X  +  X,  )2 

n  4R.  v  rn  en  tn'  v  an  tn' 

m  t—  — J 


(5) 


and 


W  = 


‘tn 

2w  . 

n  n 


Before  proceeding  with  the  applicability  of  equation  5  to  array  antennas,  it  will 
be  helpful  to  establish  a  physical  understanding  of  some  of  the  terms  involved. 

til 

a.  i  is  the  current  in  the  n  element.  This  current  is  a  function  of  the  anten- 
n 

na  element  radiation  pattern.  Further,  the  pattern  used  to  modify  i  (ex¬ 
pressed  as  g  (u)),  must  be  considered  the  pattern  of  the  element  as  it  exists 

in  the  array;  i.e. ,  all  mutual  coupling  effects  in  pattern  modification  are 
included. 

b.  If  the  array  is  an  amplitude  tapered  array,  the  values  of  i  must  be  multi¬ 
plied  by  the  tapering  coefficients,  aR.  If  the  network  provides  tapering  in 
a  lossless  manner,  this  must  be  accounted  for  in  the  basic  expression. 

c.  Ztn  =  Rtn  +  jXtn  is,  in  general,  not  a  function  of  scan  in  steerable  arrays. 

The  load  impedance  is  usually  adjusted  to  provide  the  best  average  match 
over  the  e.  ctor  of  scan. 

d.  In  a  large  array  of  equally  spaced  (or  near  equally  spaced)  identical  elements, 
the  element  pattern  function  gg(u)  and  the  element  impedance  function  Z&n  can 

be  considered  a  constant  for  all  elements.  The  large  array  stipulation  implies 
that  the  edge  effects  are  negligible. 

e.  The  power  available  in  the  receiver  (W)  is  always  a  maximum  when  the  re¬ 
ceiver  impedance  (Z.)  is  the  complex  conjugate  of  the  antenna  impedance  (Z  ). 

t  ci 

The  current  induced  by  the  incident  wave,  however,  is  a  function  of  the  termi¬ 
nal  impedance  of  the  array  element  which  is,  in  turn,  a  function  of  the  self  and 
mutual  impedances  in  the  array.  It  has  been  shown,  Referei.je  (8)  ,that  the  input 
impedance  of  the  n^  element  is  related  to  the  self  and  mutual  impedances  by 
the  following  relationship 

m 

Z  =  Z  +  y  Z  EXP  (0  )  (6) 

an  nn  ^  ^  np  '  p'  '  ' 
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vn  -  VOLTAGE  INDUCED  IN  THE 
n01  ELEMENT 

ln  -  CURRENT  IN  nth  ELEMENT 

Zan  -  INPUT  IMPEDANCE  OF  THE 
11th  ELEMENT 

Zta  -  THE  COMPONENT  OF  Zt 
TERMINATING  THE  11th 
ELEMENT 


Array  Equivalent  Circuit 


til  til 

where  Z  is  the  mutual  impedance  between  the  n  and  the  p  element  and  0 

th  th  ^ 

is  the  phase  of  the  current  in  the  p  n  element  relative  to  that  in  the  nn element. 

It  must  be  noted  that  Z  is  a  function  of  antenna  element  spacing  and  therefore 

must  be  considered  a  variable  whenever  element  spacing  is  employed  as  a  means  of 
controlling  the  array  radiation  pattern. 

Utilizing  these  comments,  it  is  now  possible  to  establish  a  generalized  expression 
of  the  signal  power  which  is  available  in  the  receiver. 

From  e. ,  above 

(W  )  =  i2  R  (7) 

'  n  max  n  rn  '  ' 

From  a. ,  b. ,  and  e. ,  above, 

g  (u)  a 

i  =  -^2 - -  EXP  <0n)  (8) 

an 

th 

where:  g  (u)  is  the  element  pattern  function  of  the  n  element,  a  is  the  relative 
en  n 

voltage  amplitude  of  the  n  element  in  an  amplitude  tapered  array,  Z  is 

the  input  impedance  of  the  n®1  element  given  by  equation  6.  an 

th 

0n  is  the  relative  phase  of  the  nn  element. 

The  maximum  value  of  signal  in  the  n^  receiver  is  therefore 


(W  ) 

v  n'max 


SWu>f  an  Rr„ 


EXP  (0n) 


or  the  total  receiver  power  becomes 


=  y.  I  Se„<“>  1 2  1  an  EXP  I  2  R 

* — t - x -  vr 


the  afi  term  in  the  denominator  is  required  to  meet  the  lossless  (efficient)  tapering 
condition  of  comment  b . 

The  signal  power  expression,  equation  9,  is  general  in  that  the  power  is  a  func¬ 
tion  of  u  for  any  position  of  the  array  beam  specified  by  0r.  In  a  large  fixed  beam 

array  the  conditions  of  comment  d  are  met  and  the  signal  is  a  maximum  when  the 
array  beam  is  aligned  with  the  direction  of  the  source.  In  the  fixed  beam  broadside 
array  the  signal  power  becomes 


W(Uq)  = 


IWI^t  I?  an  I 


I 

n  n 
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Allen,  Reference  (9),  has  shown  that  the  si,'"vnation  of  equation  10  is  simply  the 
total  number  of  array  elements  (N)  multiplied  lp,  an  amplitude  taper  efficiency  factor 
(rj)  so  that 


W(u0) 


Rt^N 


(11) 


Equation  11  states  that  the  total  signal  power  present  in  the  receiver  is  the  signal 
power  per  element  multiplied  by  the  total  number  of  elements  in  the  array. 

2.  NOISE  POWER 

In  the  black  box  concept  depicted  in  Figure  9 ,  there  are  three  conceivable 
sources  of  noise  power  which  can  be  observed  at  the  output  terminal  (dissipated  in  the 
load  resistor  Rt).  These  three  sources  have  been  identified  by  Waldman  and  Wooley, 
Reference!  Q  and  are  classified  as  follows: 

a.  External  Noise 

Noise  which  is  external  to  the  system  but  is  transferred  to  the  terminals 
as  power  collected  by  the  array.  This  power  is  a  function  of  antenna  gain 
(G  )  and  the  spatial  distribution  of  noise  power  density  ( rj ).  This  is  correlated 

noise  power. 

b.  Antenna  Noise 

Noise  generated  in  the  antenna  is  taken  to  be  the  noise  generated  in  each  an¬ 
tenna  element  and  transferred  to  the  black  box  output  by  means  of  the  network. 
This  is  uncorrelated  noise  power. 

c.  Black  Box  Noise 

The  internal  noise  power  of  the  network  is  transferred  to  the  output  and  can  be 
evaluated  as  an  equivalent  circuit  noise  factor. 

The  total  noise  power  (N^)  is  considered  to  be  the  sum  of  these  incremental  contri¬ 
butions.  Each  contribution  must  be  considered  separately  in  terms  of  the  influencing 
factors. 


A  generalized  analysis  of  noise  in  active  and  passive  circuits  is  contained  in 
Appendix  B.  An  expression  of  the  noise  power  contributed  by  a  network  is  given  by: 


AN 

where:  K 

T 

e 

B 

G 


KTgBG  equation  7  in  Appendix  B 
Boltzmann's  Constant 

Effective  temperature  of  the  network  referred  to  the  input  terminals 

Bandwidth 
Network  Gain 


The  noise  contributions  of  the  three  basic  sources  can  be  expressed  in  terms  of  this 
expression. 
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2.1  EXTERNAL  NOISE 


The  component  of  noise  introduced  to  the  system  through  the  antenna  terminals 
is  a  function  of  antenna  gain  (G  )  and  the  noise  power  density  in  space  (77) .  Since  this 

cl 

noise  is  correlated,  it  is  also  dependent  on  the  gain  (or  loss)  of  the  black  box  when 
referred  to  the  output  terminals. 


N 


x 


where  T 

x 


KBTxGr 

mja  Ga(u)r?(u)dfi 


(12) 


G 


r 


G 


r 


Effective  gain  of  the  black  box  network 


z 


G 

n 


(G^  is  the  incremental  gain  of  the  n  "  element  transferred  to  the  output 
terminal.) 


Since  the  antenna  gain  is  proportional  to  the  received  signal  power,  the  gain  of 
a  fixed  beam  broadside  array,  large  enough  to  avoid  edge  effects,  is  of  the  form 


G 


a 


R. 


(13) 


and  the  noise  power  density  is  related  to  space  temperature  by 
r?  (u)  =  Ts(u)  KB 


(14) 


Combining  equations  12,  13,  and  14,  the  external  noise  power  present  in  the  receiv¬ 
er  is  given  by 


N  = 
x 


47T  kB  R.  G 
t  r 


Wu) 


TS<U> 


z 

n 


£  a  ej7rnU 
n  n 


d^ 


(15) 


The  space  temperature  term  [T  (u)]  is  a  function  of  the  antenna  orientation  and 

s 

frequency.  The  magnitude  of  this  term  is  affected  by  ground  or  sea  temperature, 
noisy  stars,  the  sun,  the  moon,  atmospheric  absorption,  and  general  sky  or  galatic 
noise . 


The  range  of  space  temperature  is  shown  in  Figure  12  The  two  solid  curves 
represent  the  maximum  and  minimum  noise  temperature  which  will  be  present  in 
the  antenna.  The  shape  of  the  curves  indicate  the  "window"  in  the  electro-magnetic 
spectrum(l  to  10  GHz)  between  cosmic  noise  and  molecular  absorption. 

In  order  to  completely  specify  the  external  noise  conditions,  the  gain  and  noise 
temperature  conditions  must  be  known  for  the  sector  of  space  in  which  the  system  is 
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SPACE  TEMPERATURE  (K°) 


Note:  Curves  based  on  data  of  J.  C.  Green  and  M.  T.  Lehenbaum, 
Letter  in  Microwave  Journal .  Vol.  2,  pp.  .13,  14 
October  1959. 


F  igure  1  d 


Noise  Temperature  vs.  Frequencv 


to  operate.  Since  this  is  likely  to  vary  within  a  particular  application,  maximum 
and  minimum  limits  can  he  established  as  an  estimate  oi  performance. 


2.2  ANTENNA  NOISE  POWER  (N  ) 

Each  antenna  element  is  a  source  of  noise  with  a  related  noise  temperature  T0n. 
The  total  effective  noise  temperature  of  the  antenna  is  given  by 


where  G  is 
n 


T 
'  a 


2t 

n  en 


G 
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r 

the  gain  of  the  network  from  the  n 


th 


terminal  to  the  output 


N 
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KBT  G 
a  r 


v 

=  KB  2-j  T_„G 


en  n 


(16) 


2.3  NETWORK  NOISE  POWER  (N  ) 

r 

The  effective  noise  temperature  of  the  black  box  network  is  T  ,  therefore  the 
internal  noise  present  at  the  output  terminal  is 

N  =  KBT  G  (17) 

r  r  r 

The  total  noise  power  present  at  the  output  terminal  is  simply  the  sum  of  the 
three  contributions  and  is  given  by 


or  N,  =  KBG  (T  +  T  +T  ). 

t  r  x  a  r 


(19) 


Before  taking  the  ratio  of  equations  10  and  18  to  determine  the  signal-to- noise 
conditions  in  the  array,  some  observations  covering  the  noise  expression  in  equa¬ 
tions  18  and  19  are  in  order.  The  stated  purpose  of  this  study  is  to  determine  the 
effects  of  element  spacing  and  weighting  on  the  signal-to-noise  ratio  of  the  system. 
It  is  noted  that  if  the  array  is  comprised  of  a  fixed  number  of  elements,  the  only 
term  in  equation  19  which  is  affected  by  weighting  and  spacing  is  the  first  term  T  . 

If  the  number  of  elements  is  not  fixed,  the  Ta  term  may  be  affected  but  if  the  array 
elements  are  all  assumed  to  be  at  the  same  temperature  (Te),  then  Ta  =  Tg  since 
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A  logical  approach  in  improving  the  array  signal-to-noise  ratio  is  to  consider  the 
ratio  of  output  signal  to  output  noise  resulting  from  the  external  noise  power  density 
(?]).  Assuming  the  array  beam  is  oriented  to  receive  maximum  signal,  the  signal-to- 
noise  ratio  which  is  affected  by  element  spacing  and  weighting  is  given  by  taking  the 
ratio  of  equations  10  and  15,  keeping  in  mind  that  the  signal  at  the  output  terminal  is 
WGr>  After  cancellation  of  common  terms,  this  signal-to-noise  ratio  becomes 


(20) 


The  significance  of  equation  20  can  be  generally  stated  in  words.  The  numerator 
represents  the  maximum  value  of  the  antenna  element  gain  function  multiplied  by  the 
maximum  value  of  the  array  function  which  is  equivalent  to  the  element  gain  times  the 
number  of  elements  in  the  array.  The  denominator  is  seen  to  be  a  constant  times  the 
intergal  (or  summation)  of  the  element  gain  function  times  the  noise  distribution  times 
the  array  pattern  function .  It  is  this  array  pattern  function  which  has  been  analyzed 
in  Section  III. 

From  equation  20,  it  can  be  seen  that  the  signal-to-noise  is  dependent  on  a  num¬ 
ber  of  factors,  including  the  following: 

a.  Array  gain 

b.  Noise  temperature  distribution 

c.  Antenna  element  pattern  (as  observed  in  the  array,  including  coupling  effects) 

d.  Array  pattern 

e.  Beam  pointing  angle  (uq)  . 

In  order  to  evaluate  any  array  for  S/N,  these  parameters  must  be  known  or  assumed. 

3.  MAXIMIZATION  OF  SIGNAL-TO-NOISE  RATIO 

Considerable  effort  in  the  past  has  been  devoted  to  the  process  of  increasing  the 
signal-to-noise  ratio  in  array  antenna  systems.  Neglecting  the  effect  of  external 
noise;  i.e. ,  assuming  a  uniform  environment,  these  efforts  have  been  directed  toward 
(1)  the  minimization  of  uncorrelated  noise  and  (2)  the  maximization  of  correlated  signal. 

The  former  approach  has  been  pursued  by  Linder,  Reference (11), and  expanded  by 
Cheng  and  Tseng,  Reference ( 12) .  The  process  requires  the  representation  of  the  array 
cross-correlation  coefficients  in  a  matrix  which  can  be  minimized  by  the  application  of 
Lagrange  multipliers  provided  the  quantity  of  array  elements  is  small.  This  is  equiva¬ 
lent  to  an  inversion  of  the  matrix  which  yields  a  solution  in  element  weighting  or  ampli¬ 
fier  gain.  The  results  obtained  by  Linder  indicate  that  a  significant  improvement  in 
signal-to-noise  ratio  can  be  achieved  with  closely  spaced  elements  (less  than  one-half 
wavelength)  provided  the  severe  amplitude  weighting  can  be  maintained. 

The  maximization  of  correlated  signed  is  an  attempt  to  increase  the  directive  gain 
of  the  antenna  array.  This  technique  has  been  investigated  extensively,  References  (13)  thru  (18). 
In  a  linear  array,  it  has  been  shown  by  Tai,  Reference  (19),  that  the  directivity  can  be 
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increased  by  close  element  spacing  and  a  severe  amplitude  taper.  The  results  of  Tai 
and  the  mathematical  formulation  exactly  parallels  that  of  Linder  even  through  the 
approach  to  the  problem  was  based  on  different  initial  conditions. 

The  increase  in  array  directivity  or  signal-to-noise  ratio  associated  with  close 
element  spacing  and  severe  amplitude  tapering  is  commonly  referred  to  as  the 
"super-gain  effect".  Lo,  et.  al.,  Reference  (20).  relates  the  super -gain  effect  to  a 

quality  factor  Q.  For  any  set  of  conditions  there  is  practical  limit  on  the  Q-factor 
which  determines  the  degree  of  super-gain  permissible.  The  Q-factor  is  essentially 
an  expression  of  the  bandwidth  sensitivity  and  efficiency  of  the  super-gain  array. 

The  super-gain  effect  is  illustrated  in  Figures  13  and  14  .  The  curve  in  Figure  13 
depicts  the  variation  in  array  directivity,  associated  with  a  10-element  array, 
as  a  function  of  element  spacing.  It  is  noticed  that  for  element  spacing  less  than 
0.5  wavelengths  it  is  possible  to  achieve  the  super-gain  condition  through  applica¬ 
tion  of  the  prescribed  optimum  illumination  function.  The  same  information  is 
shown  in  Figure  14  where  the  array  length  is  fixed  at  five  wavelengths  while  the 
number  of  elements  comprising  the  array  is  allowed  to  vary. 

Figure  13  shows  that  the  highest  directivity  which  can  be  achieved  from  a  10- 
element  array  is  about  12  DB  corresponding  to  element  separation  of  about  0.9 
wavelengths.  Figure  14  shows  that  the  highest  directivity  which  can  be  obtained  in  a 
5 A  aperture  without  super-gain  is  10  DB  and  this  is  approached  with  a  minimum  of 
6  elements  also  separated  about  0.  9  wavelengths. 

A  note  concerning  the  practicality  of  the  super- gain  effect  is  in  order.  Although 
often  stated,  the  limitations  on  achieving  super-gain  cannot  be  over  emphasized.  The 
super-gain  condition  is  achieved  by  exciting  the  array  elements  with  a  severe  ampli¬ 
tude  taper  which  reverses  phase  on  an  element-to-element  basis.  This  condition  is 
extremely  difficult  to  establish  and  maintain.  The  degree  of  super-directivity 
possible  is  dependent  on  the  element  spacing.  As  the  elements  are  brought  closer 
together,  interelement  coupling  becomes  stronger,  making  adjustment  difficult.  The 
excitation  amplitudes  become  so  critical  that  an  error  of  one  percent  can  completely 
destroy  the  super-gain  effect.  In  application,  therefore,  super-directivity  is  limited 
to  the  narrow  band  situation  where  precise  control  is  available  to  achieve  a  decrease 
in  array  beam  width. 

Although  the  terms  super-gain  and  super-directivity  are  used  interchangeably, 
the  latter  is  correct.  As  the  directivity  is  increased  beyond  that  obtainable  with 
uniform  illumination,  the  efficiency  is  degraded  accordingly,  due  to  high  currents 
(or  gain  requirements  in  amplification)  required  to  produce  the  optimum  illumina¬ 
tion  function. 

Although  the  super-directive  technique  is  not  feasible  as  a  means  of  achieving 
high  gain  from  small  apertures,  some  improvement  in  signal-to-noise  ratio  can  be 
obtained  if  other  aspects  of  system  performance  can  be  sacrificed.  For  example, 
low  noise,  variable-gain  amplifiers  might  be  used  in  a  receiving  array  to  produce 
an  improved  signal-to-noise  ratio  through  gain  adjustment.  Improvement  of  this 
nature  appears  to  be  limited  to  about  10  percent  under  realistic  conditions. 

The  study  of  signal-to-noise  ratio  by  Cheng  and  Tseng,  Reference  (12),  has 
provided  a  general  formulation  of  the  optimization  process.  It  is  noted  that  there 
are  eight  basic  factors  which  affect  the  ultimate  signal-to-noise  ratio  of  a  fixed 
array  antenna  system . 
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a.  Antenna  array  configuration  and  element  spacing 

b.  Element  radiation,  impedance,  and  frequency  characteristics 
(including  coupling  effects) 

c .  Source  direction  relative  to  the  array 

d.  Receiver  or  filter  response 

e.  Signal  spectrum 

f  .  Noise  distribution  in  space  and  frequency 

g.  Internal  noise  considerations 

h.  Element  weighting  or  amplification  factors. 


It  has  also  been  shown  by  Lo,  et.al. ,  Reference  (20),  that  if  the  noise  distribu- 

in  space  and  frequency  is  uniform,  the  conditions  yielding  maximum  signal-to- noise 
ratio  and  maximum  directivity  ire  identical. 


The  directivity  of  any  antenna  is  defined  at  the  ratio  of  the  peak  radiated  power 
per  unit  solid  angle  to  the  average  radiated  power  per  unit  solid  angle  where  the 
averaging  process  includes  all  space  angles.  Mathematically  stated,  using  spherical 
coordinates,  this  is 
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P  sin©  dG 
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P  (0,  <p)  sin©  dG  d$ 


(21) 


which,  since  P  is  constant,  equation  21  is 
max 


4?r  P 
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(22) 


For  a  broadside  linear  array  of  N  discrete  elements  with  feed  coefficients,  A  , 
aligned  on  the  z-axis, 


P(G,0) 

where 

E  (©  ,<p) 
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(23) 

(24) 
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If  the  element  pattern  function  is  taken  as  unity  so  that  G(Q,<p)~  1,  the  linear 
array  consists  of  N  isotropic  radiators.  Further,  if  the  feed  coefficients  are 
cophasal,  the  beam  maximum  will  occur  at  9-tt/2.  The  radiation  pattern  can 
then  be  normalized  to  unity  so  that  P  =1  and  equation  (24)  becomes, 
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(  25  ) 


Using  equation  (  25  )  with  equation  (  22)  gives, 
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sin  0d9 


Maximizing  the  directivity  amounts  to  minimizing  the  integral 
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An  alternate  way  of  writing  the  factor  which  is  the  square  of  the  absolute 
value  in  equation  (  27 )  is, 
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Then, 
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By  letting  x  =  cos  9,  the  integral  is  simplified  to 

r1 


(n-ro) 


XdX 


36 


Or,  upon  performing  the  indicated  integration, 
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(  29) 


In  the  general  case  the  value  of  this  term  is  dependent  on  the  antenna  ele¬ 
ment  gain  function  ( G(Q,cp )  in  equation  24),  the  element  amplitudes  (A  ),  the  an¬ 
tenna  element  spacing  (d),  the  number  of  elements  in  the  array  (N)  an$  the  operating 
frequency  (f  =  cA). 

The  expression  "I”  is  defined  for  a  uniform  array  of  isotropic  sources.  The 

equation  can  be  made  general  by  replacing  the  sin  x/x  term  with  the  generalized 

cross  correlation  coefficient  between  the  n^1  and  m”1  elements  (R  ). 
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(30) 


It  is  clear  that  R  is  dependent  on  the  distance  between  the  two  elements  as 
n,  m 

well  as  the  applied  frequency  conditions. 


The  condition  for  maximum  directivity  is  found  by  differentiating  equation 
(30)  with  respect  to  each  of  the  amplitude  c  oefficients  and  setting  the  result  equal 
to  zero. 


or 
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Equation  31  represents  a  set  of  N  simultaneous  equations  where  k  = 
1,2,3,  — ,  N.  Setting  each  of  these  equations  to  zero  and  solving  for  yields 

the  amplitude  coefficients  required  to  maximize  the  array  directivity. 

A  matrix  inversion  method  can  be  used  to  solve  for  the  amplitude  co¬ 
efficients  of  equation  31  .  This  expression  can  be  rewritten  in  the  form 
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which  can  also  be  expressed  in  matrix  form 

[r]n,n  Hn,i  =0n,i 

where  [k]  is  a  column  vector  of  identical  elements  equal  to  the  right  hand  side  of 
equation  32  .  Premultiplying  by  the  inverse  matrix  of  [Rj  gives 

Mn.i  =  [d1n,n  Hn.i  'Vhere[D]  =  M 

The  amplitude  coefficients  are  then  proportional  to  the  sum  of  the  matrix  elements 


Ak  a  £  Dk,j  K  =  1,2,3,  , N 


In  a  large  array  the  matrix  inversion  required  for  solution  becomes 
difficult.  In  a  symmetrical  array,  the  matrix  inversion  can  be  simplified  by 
shifting  the  reference  from  the  centerline  of  the  array  to  the  center  of  one  side. 
This  is  accomplished  by  redefining  the  correlation  coefficient 


)3  -  R  +  R  XTj, 

n,  m  n,m  n,  N+l  -  m 


This  transformation  holds  for  both  odd  and  even  arrays.  Making  this  substitution 
in  equation  31  the  new  form  is 
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The  solution  to  equation  33  is  identical  with  that  described  above  for  equation  31 
The  advantage  is  that  the  new  matrix  [P  ]  is  one  half  the  order  of  the  matrix  [R] 
which  simplifies  the  inversion  required  for  solution. 


For  very  large  arrays  (  N  >  100  )  even  the  simplified  form  of  equation 
31  is  difficult  to  solve  by  means  of  matrix  inversion.  For  this  reason  an  iterative 
technique  has  been  developed  which  is  described  in  Section  VII.  Four  forms  of  the 
basic  equation  expression  31  are  used  to  implement  the  iterative  process 
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(  37) 


Equations  34  and  35  are  based  directly  on  equation  31  while  36  and  37 
are  derived  from  its  counterpart,  equation  33  .  The  equation  to  be  used 
depends  on  the  nature  of  the  ar  ray.  Equations  36  and  37  apply  only  to  symmetrical 
arrays  but,  in  general,  produce  a  higher  rate  of  convergence. 
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The  iterative  process  begins  by  assuming  values  for  the  excitation 

coefficients  such  as  A  =1  for  all  n.  After  the  A,  terms  have  been  computed, 

n.  k 

these  terms  are  used  to  replace  the  A^  terms  and  the  process  is  repeated  until 

the  change  is  negligible.  The  rate,  of  convergence  is  somewhat  dependent  on 
which  elements  are  computed  first. 

Even  though  the  iterative  process  does  converge,  there  is  no  assurance 
that  the  values  obtained  are  optimum.  It  is  usually  desirable  to  use  more  than  one 
of  the  above  expressions  as  a  final  check  on  the  solution. 


4.  MAXIMIZATION  OF  ARRAY  DIRECTIVITY  THROUGH  ELEMENT  SPACING 


The  process  of  optimizing  array  directivity  through  control  of  the 
antenna  element  location  is  directly  analogous  to  the  process  described  above 
where  the  amplitude  of  excitation  was  used  as  the  controlling  factor.  The  basic 
equation  30  is  again  differentiated,  but  this  time  with  respect  to  d^  rather  than 

Ak‘ 

Before  differentiating  it  is  convenient  to  rewrite  equation  30  for  symmetric 
array  with  the  center  of  the  array  used  as  a  reference.  Under  these  conditions 
the  denominator  of  equation  30  becomes 
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for  an  even  array,  the  numerator  of  equation  30  is 
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For  an  array  having  an 'I  odd  number  oi  elements,  the  double  sum  is 

P  P  P 


A  +  2A 
o  o 


2,  A  (R  +  R 

n'  o,n  0, -n 


)  +  2  Z  Z  A  A 

'  n 


n=l 


n=l  m=l 


n  m 


R  +  R 
n,m 


n,  -nil 
(39) 


or 


A  +  1A  Z,  A  R  +2 


n=l 


n  o,n 


Z  Z  A  A  Tr  +R 

,  ,  n  m  I  n,m  r.,-m 

n=l  m=l  L 


40 


It  is  noted  that  the  summations  for  the  odd  array  reduce  to  the  even  array  formula¬ 
tion  by  simply  letting  Aq  equal  rero.  Thus,  a  solution  of  the  odd  array  case  can  be 

applied  to  an  array  having  an  even  number  of  elements  quite  simply.  The  derivation 
is  performed  for  the  odd  array  for  this  reason. 
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Substituting  equations  38  and  39  into  equation  30  for  the  odd  element  con- 
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In  this  form  differentiation  with  respect  to  dx  is  straightforward. 
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This  expression  can  be  simplified  considerably.  Let  R'(d)  be  the  first  derivative 
of  the  correlation  coefficient  with  respect  to  the  distance  between  two  elements. 
Then,  since  R'(d)  is  an  odd  function, 

R'(-d)  -  -R'(d), 

and  it  becomes  apparent  that  the  crder  of  the  subscripts  is  significant.  The  con¬ 
vention  adopted  here  is 


R’(d)  =  R'(dn-dk) 
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Using  this  convention  and  notation,  the  derivatives  indicated  in  equation  41  are 
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Applying  these  simplifications,  equation  41  becomes 


n=l 


where  the  n  ^  k  restriction  of  equation  41  has  been  removed.  Equating  this 
derivative  to  zero  and  dividing  through  by  the  constants  gives  the  condition  of 
maximizing  array  directivity  through  element  position  control. 
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The  appearance  of  equation  43  is  less  formidable  if  a  new  variable  is 
defined.  Let 


*k.„  -  -R,<W  - ey-V 

Then  equation  43  becomes 
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n=l 

where  k  =  1,2, ....  P,  since  there  are  P  possible  element  pair  locations.  Fixing 

some  of  these  element  pairs  amounts  to  imposing  a  constraint  on  some  of  the  pair 

positions.  Where  d,  is  constrained  to  a  specified  fixed  position,  y,  =0  and 
k  k,  n 

R\  =  0.  Thus,  k  =  1,2, . .  .  ,  L,  where  L  is  the  number  of  degrees  of  freedom 

K)  0 

in  the  system  of  equations.  A  good  example  of  the  practicality  of  fixing  an  element 
pair  position  is  to  specify  the  position  of  the  end  e^ment  pair,  thereby  limiting 
the  overall  aperture.  Equation  45  written  in  matrix  form  is 

Wl,p  Hp>1  -  Ao  [r,]Li1  (46) 


(  44) 


(45) 


Although  the  form  of  equation  46  is  identical  with  the  amplitude  weighting 
solution,  the  problem  of  determining  the  actual  element  pair  positions  from  the 
equation  presents  some  difficulty.  This  arrises  from  the  fact  that  a  solution  is 
not  generally  unique.  This  non-ur.iqueness  may  result  in  many  minor  maximum 
values  of  array  directivity.  Since  the  absolute  maximum  directivity  is  normally 
desired,  it  is  apparent  that  equation  46  is  a  necessary  but  not  a  sufficient  condition. 
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The  solution  to  equation  43  might  be  identified  as  an  experimental  process 
similar  to  one  which  could  be  encountered  on  an  antenna  test  range.  The  first  step 
in  experimental  solution  would  be  to  determine  the  spacing  required  to  obtain  the 
maximum  directivity  fcr  an  equally  spaced  array.  However,  the  maximum  directivity 
for  an  array  of  unequally  spaced  elements  must  be  higher  than  that  of  the  optimum 
array  using  equal  element  spacing.  Thus,  the  equally  spaced  array  serves  as  a 
starting  point.  The  next  step  is  an  attempt  to  increase  the  directivity  by  adjusting 
the  positions  of  each  pair  of  elements,  one  pair  at  a  time.  The  process  can  be 
repeated  until  there  is  no  perceptable  movement  of  the  element  positions.  This  pro¬ 
cedure  can  be  applied  mathematically  using  the  appropriate  equations. 


For  a  given  set  of  amplitude  coefficients,  the  array  directivity  can  be  de¬ 
termined  for  equal  element  spacing  by  applying  equation  26  .  Then  the  necessary 
condition,  represented  by  equation  43  ,  can  be  applied.  For  computation  purposes, 
equation  43  can  be  written  in  the  form, 
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where  the  solution  for  d^  is  such  that  F(d^)  =  0.  To  obtain  this  solution,  a  modifica¬ 
tion  of  Newton's  method  for  finding  roots  can  be  applied.  The  solution,  according  to 
Newton's  method  is 


dk  -  dk  -  F«y  /  F,«y 


(48) 


Since  the  correction  on  -  F(d^) ,  can  be  quite  large,  some  upper  limit 
on  the  correction  should  be  imposed.  The  derivative  of  equation  47  is 

P 

F'<dk>  '  Vk,o  *  2AkRk,-k  +  £  An<R' k.n  +  Rk,-n>  (49) 
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After  using  the  modification  of  Newton's  method  for  finding  the  k  element 
pair  locations,  the(k±l )-th  pair  can  be  determined.  By  going  through  the  array 
several  times  in  this  manner,  any  possible  improvement  in  the  array  directivity 
can  be  found.  A  difficulty  with  this  approach  is  that  there  is  no  guarantee  that  the 
directivity  will  be  the  absolute  maximum,  although  it  will  be  higher  than  the  equally 
spaced  array. 


As  an  example  of  the  results  which  can  be  obtained  through  element  control, 
two  array  conditions  are  presented  in  Table  3  .  The  relative  element  spacing  is 
compared  with  that  of  a  comparable  uniform  array,  and  the  directivity  values  are 
presented  indicating  the  degree  of  improvement  available. 
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TABLE  3.  RESULTS  OF  ELEMENT  SPACING  COMPUTATIONS 


6  Elements  —  Uniform  Amplitudes 

Element 

Element  Spacing  for 

Number 

single  frequency 

UNIFORM 

OPTIMUM  DIRECTIVITY 

1 

157.8 

161.7 

2 

473.4 

482.3 

3 

788.9 

778.8 

S/N 

9.788  DB 

9.844  DB 

20  Elements  —  Uniform  Amplitudes 

Element 

Element  Spacing  for 

Number 

single  frequency 

UNIFORM 

OPTIMUM  DIRECTIVITY 

1 

165.0 

173.4 

2 

495.0 

520.0 

3 

824.9 

866.2 

4 

1154.9 

1212.0 

5 

1484.9 

1556.9 

6 

1814.9 

1900.4 

7 

2144.9 

2241.9 

8 

2474. 8 

2580.0 

9 

2804.8 

2911.8 

10 

3134.8 

3217.2 

S/N 

15.410  DB 

— 

15.623  DB 
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SECTION  V 


WIDEBAND  EFFECTS 


The  array  design  considerations  and  the  signal,  noise,  and  directivity  analysis 
presented  in  Sections  III  and  IV  were  limited  to  the  single  frequency  condition.  Most 
systems  do  not  operate  at  precisely  a  single  frequency,  i.e.  with  a  continuous  wave 
input,  but  rather  operate  over  a  range  of  frequencies  defined  by  input  waveform .  noise 
conditions,  tunability  of  equipment  or  '>ther  component  characteristics.  Unless  the 
bandwidth  is  very  wide,  the  single  frequency  assumption  provides  an  acceptable  base 
from  which  to  determine  approximate  system  performance. 

In  this  section,  the  narrow  band  assumption  is  removed  and  the  effect  of  wide¬ 
band  signals  on  array  performance  is  determined  in  terms  of  aperture  effects  and 
signal -to -noise  ratio  effects.  A  possibility  of  utilizing  the  wideband  condition  as  a 
means  of  improving  array  performance  is  also  presented. 

1.  APERTURE  EFFECTS 

The  general  expression  for  the  radiation  pattern  of  an  array  of  N 

elemeutc  is  given  by 

N  f 

E(u,t)=  e^2,rft  ^  Ane^nr>u  (50) 

n  -  1 


where 


f  =  the  frequency  of  the  incoming  wave 
fQ  -  an  arbitrary  "center"  frequency 

u  =  |^  sin  0,  d  is  the  element  spacing 
-  the  wavelength  of  the  center  frequency 

9  =  the  observation  angle  relative  to  boresight 
Afi  =  the  array  element  amplitude  coefficients 

In  the  analysis  of  narrow  band  arrays,  it  is  often  convenient  to  neglect  the  term  with 

a  statement  that  it  is  "understood"  to  be  a  multiplier  of  all  expressions.  This  is  valid 
only  if  the  applied  signal  is  independent  of  time,  i.e.  a  CW  signal.  If  more  than  one 
frequency  is  applied  to  the  array  simultaneously,  a  component  of  E(u,t)  will  result  for 
each  frequency,  so  that  equation  50  becomes 
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where  is  the  amplitude  weighting  of  the  k**1  frequency.  If  the  frequencies  are  related 
so  that 

fk  =  fo+kAf  k  =  0,  ±1,  ±2,  ±3,  — ,  ±  (52;. 

(Af  is  a  small  difference  frequency;,  equation  52  can  then  be  rewritten 


E(utt)  =  f  FkAnei **■£-■,*  (53) 

k  =1  n  -  i 

Equation  53  represents  a  repeating  waveform  with  a  period  T  =  ^  .  The 

spectral  lines  comprising  the  waveform  are  characterized  by  the  phase  and  amplitude 
of  the  complex  form  ^  ej2ff(fo  <  kAf)tt 

k 


In  application,  it  is  common  to  separate  signal  information  from  its  RF  car¬ 
rier  by  means  of  detection.  The  detection  circuit  provides  the  means  by  which  the 

e^^o*  term  in  equation  53  is  removed.  Since  this  carrier  is  a  time  varying  signal, 
the  time  constant  of  the  detector  circuit  is  extremely  important  The  requirements  for 
good  detection  are:  the  detector  time  constant  must  be  long,  as  compared  with  the  RF 
cycle,  and  short  as  compared  with  the  desired  video  waveform.  These  two  require¬ 
ments  may  not  be  compatible  where  extremely  large  bandwidth?  are  involved  but,  for 
purposes  of  analysis,  perfect  detection  is  assumed. 


Removing  the  RF  component  from  equation  53  yields  the  detected  output. 


E(u,  t) 


j7Tk(2Aft  +  npu)| 
k  o  J 


(54) 


Equation  54  is  based  on  the  excitation  of  an  array  containing  a  fixed  number 
of  elements  N,  with  a  fixed  number  of  spectral  components  K.  It  is  possible  to  establish 
a  comparable  expression  wherein  the  illumination  function  and  the  frequency  spectrum 
are  continuous.  The  video  response  corresponding  to  the  continuous  function  is 

fo  +  6  1 

E<z,t>  -  /  \  "  F(0  G(X)  ei2"i-  **  e)2»(t  -  fQ)t  dxdJ  <55) 

o  z 


where 


sin  0,  A  -  the  aperture  in  wavelengths  (center  frequency) 
2  6  =  the  bandwidth 

G(x)  =  the  aperture  amplitude  distribution 
F(f)  =  the  frequency  amplitude  distribution 


In  the  manner  of  equations  54  and  55  the  combination  of  continuous  and  dis¬ 
crete  functions  can  be  expressed. 
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Discrete  aperture  with  continuous  spectrum 

f  +  6 

I  F  if)  ei2'  <f  -(o> 1 


E(u,t)  = 


f  -6 
o 


Continuous  aperture  with  discrete  spectrum 

1 

f~  2  j27T  # —  ZX 

E(z,tj  =  J ^  G(x}  e3  iQ 

~  z 


N 

£ 
n  =  I 


Ane37mu  df 


(56) 
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k  =  1 


ux 
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Each  of  the  above  expressions  represents  a  radiation  function  which  is  variable 
in  both  space  and  time.  This  means  that  a  true  representation  of  the  wideband  array  can 
only  be  observed  in  the  time-space  coordinate  system.  An  example  of  the  power  re¬ 
sponse  observed  in  the  time-space  domain  for  a  uniform  discrete  array,  illuminated 
with  a  uniform  discrete  frequency  function,  is  shown  in  Figure  15. 

The  three  dimensional  model  of  Figure  15  i;.  used  to  illustrate  the  distribu¬ 
tion  of  energy  as  specified  by  equation  54  .  There  are  several  characteristics  of  this 
distribution  which  are  significant: 

a.  The  principal  plane  pattern  along  the  time  axis  (u  =  0),  follows  the 

x 

variation,  but  along  the  angular  axis  (tf  =  0)  decays  rapidly 
to  an  insignificant  level. 

b.  A  peak  off-axis  response  occurs  along  straight  lines  displaced  in  both 
time  and  space. 

c.  The  far  out  sidelobe  response  along  these  displaced  lines  is  6  DB  below 
the  peak  sidelobe  response  of  the  same  array  excited  with  a  CW  signal. 

An  explanation  of  these  observations  can  be  gleaned  from  a  close  examination  of 
an.  ideal  time  waveform  defined  by  the  Fourier  integral 


E.  (t)  = 


1_ 

26 


A(f) 


i2wft 

e' 


df 


(58) 


where: 

26  =  the  frequency  range 

A(f)  =  the  spectral  density 

defined  so  that  a  peak  value  is  obtained  at  t  =  t  . 

In  the  wideband  antenna  array,  the  ideal  waveform  is  modified  by  the  antenna  function 
G(f)  • 


E(t)  = 


1_ 

26 


(59) 
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Figure  15.,  Time-Space  Response  of  Wideband  Array 


The  antenna  array  sidelobe  behavior  is  defined  in  tei'ms  of  G(f) ,  but  for  many 
linear  arrays  the  peak  value  is  essentially  constant  over  a  limited  range  and  can  be  ap¬ 
proximated  by 


G(B 


f  +  f  \ 
_ I 

f0  J 


(30) 


where: 

N  =  the  number  of  array  elements 
f  =  the  center  frequency 

f  =  the  deviation  from  f 

u  -  the  far-field  space  coordinate 


Since  equation  60  represents  only  the  sidelobe  structure  of  the  array  pattern,  u  must 
be  limited  to  large  values.  Incorporating  the  cosine  identity 


cos  x 


4 


Li 


into  equation  60  and  substituting  into  equation  59  a  sum  in  the  form  of  equation 
58  is  obtained.  Utilizing  the  general  form  of  equation  58  the  result  can  be  written 


E(t) 


1 

2 


e 


i  T.  Nu 
2 


E. 

l 


TTNu  r.  Nu  \ 
2  Eiit'4f0) 


(61) 


An  understanding  of  the  sidelobe  behavior  in  the  far-out  pattern  can  be  obtained 
by  considering  the  narrow  band  and  wideband  response.  In  the  narrow  band  case  f  =  0  , 
and  the  sidelobe  pattern  is  given  by  equation  30  as  cos  ttNu/2  .  When  the  bandwidth 
becomes  large,  the  response  is  given  by  equation  61  which  consists  of  two  separate 
terms.  The  sidelobe  region  can  then  be  considered  as  the  sum  of  two  ideal  waveforms. 
In  the  far-out  region,  the  interaction  between  these  waveforms  is  very  slight  and  the 
response  contains  a  double  peak  in  time  appearing  at 

,  Nu  .  .  ,  Nd  .  „ 

t  =  t  ±  — —  or  t  -  t  ±  r —  siri  0  1 62) 

o  4f  o  2P  v  ' 

o  ° 

where: 

d  =  element  separation 
c  =■  velocity  of  light 
©  =  angular  displacement 

This  double  peak  response  is  observed  in  Figure  15  for  any  value  of  tfQ  displaced 

from  zero.  Since  the  response  is  symmetrical  about  t  =  o,  there  are  four  equal  lobes 
appearing  in  the  time-space  pattern. 

The  amplitude  of  the  sidelobe  function  in  the  region  of  slight  interaction  is 

given  by 


E(VJr)  =  1/2Ei<y  <63> 

which  occurs  at  the  peak  of  the  waveform.  The  1/2  in  equation  63  represents  the 
nominal  6-DB  reduction  observed  in  the  time-space  domain  over  that  appearing  in  the 
CW  case. 
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The  above  discussion  is  limited  to  the  behavior  of  the  array  ir.  the  time-space 
domain,  removed  from  the  area  of  the  major  lobe.  In  the  area  of  the  major  lobe,  the 
response  ij  not  so  well  defined.  The  response  of  the  uniform-array  (Figure  15)  in 
the  area  of  the  major  lobe  is  shown,  expanded,  in  Figure  16  .  It  is  noted,  in  follow¬ 
ing  the  response  along  the  u  axis,  that  the  major  lobe  broadens  ar.d  smooths  out  prior 
to  splitting  into  the  double  peak  indicated  by  equation  61  .  This  effect  is  caused  by 
severe  interaction  of  the  tv/o  terms.  The  effect  of  interaction  can  be  seen  in  Figure  17 
which  is  a  contour  diagram  of  the  response  depicted  in  Figures  15  and  16  .  The 
peak  response  is  shown  as  a  dashed  line  which  is  seen  to  meander  about  the  line  pre- 
dicteJ  by  »■  triatlor.  62  , 


sma1! 


Ig  this 


, ’slues  of  u.  The  act'.’?!  '■'’esporsst  along 
dieted  line  is  shown  in  Figure  18  ,  compared  with  peak  response  at  t  =  O.  The  dif¬ 
ference  between  these  two  curves  is  shown  to  be,  on  the  average,  about  6  DB  eorres- 
IMvJiog  ttte  ppfcdktijR.  oi  e^Mtiun  €3  .  "ffet  variation  to  ajnpl&iJ*,  as  ft  fwwlluft  (M 
u,  represents  the  degree  of  interaction  between  the  terms  of  equation  61  . 


The  time-space  response  of  an  array  is  dependent  on  the  frequency  and  ampli¬ 
tude  weighting  functions.  The  reepv,..ct  for  Is  sh^vti  in  lv.rn: 

in  Figures  17  through  28  .  The  excitation  conditions  are  indicated  in  Table  4. 


TABLE  4  .  ILLUSTRATED  TIME-SPACE  RESPONSE  FUNCTIONS 


r  LUfegram 

ill  uiiilnutk.lt 

v(f)k 

Figure 

are  No. 

Function 

Function 

17 

Uniform 

Uniform 

18 

19 

Uniform 

Cosine 

20 

21 

Cosine 

Uniform 

22 

23 

Cosine 

Cosine 

24 

25 

Uniform 

Square  Pulse 

— 

26 

Cosine 

Square  Pulse 

— 

27 

Uniform 

Trapezoidal  Pulse 

— 

28 

Cosine 

Trapezoidal  Pulse 

— 

2.  SIGNAL-TO -NOISE  RATIO  EFFECTS 


In  Section  IV  the  signal -co -noise  ratio  of  an  array  system  is  defined  in  terms 
of  a  number  of  interdependent  factors  which  represent  variables  in  the  signal  power 
and  noise  power  present  at  the  terminals  of  the  array. 


Linder,  Reference  (11)  ,  has  shown  that  the  signal-to-noise  ratio  in  a  linear 
array  can  be  expressed  as 

.  \ 

2 


S/N  = 


N 

X  n 

n  =  1 


N  N  A  a  ra 

v->  v*  AAR 

2^  n  m  n ,  m 

n  =  1  m  =  1 


(64) 


where: 

An  =  the  amplitude  of  the  nth  element 

A^  =  the  amplitude  of  the  mth  element 

R  =  the  cross-correlation  coefficient  between  the  nth 
n  rp 

’  '  and  mth  elements . 
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Figure  16„  Near  Boresight,  Response  (Expanded)  Wideband  A 


Figure  17.  Near  Boresight  Response  Contour  Plot-Uniform  Frequency 


Figure  18.  Maximum  Response  in  Time-Space  Domain 


Figure  19,  Near  Boresight  Response  Contour  Plot-Cosine  Frequency 


Figure  20.  Maximum  Response  in  Time-Space  Domain-Cosine  Frequene 


Figure  21.  Near  Boresight  Response  Contour  Plot 


Figure  22.  Maximum  Response  in  Time-Space  Domain 


Boresight  Response  Contour  Plot 


Figure  25 .  Near  Boresight  Response  Contour  Plot 


Boresight  Response  Contour  Plot 


Figure  27.  Near  Boresight  Response  Contour  Plot 


Figure  28.  Near  Boresight  Response  Contour  Plot 


In  its  general  form,  the  cross -correlation  coefficient  is  a  function  of 

a.  Separation  between  antenna  elements 

b.  Antenna  element  pattern  function 

c.  Noise  power  distribution  in  space 

d.  Signal  power  spectral  density. 

A  discussion  of  the  cross -correlation  coefficient  is  contained  in  Appendix  C. 
Graphical  representation  of  the  correlation  coefficient  is  used  to  illustrate  the  effect 
of  antenna  element  type,  frequency  bandwidth,  and  spatial  no;oe  distribution. 

The  effect  of  bandwidth  on  the  array  signal -to -noise  ratio  can  be  observed 
from  its  effect  on  the  correlation  coefficients.  In  the  linear  array  of  isotropic  source 
elements,  operating  in  a  uniform  noise  environment  at  a  single  frequency,  the  correla¬ 
tion  coefficient  is  giv^n  by 


sin  2  7T d 

R(d)  =  X  (65) 

27rd 
\ 

where  d  is  the  separation  between  elements. 

If  the  single  frequency  restriction  is  removed,  the  correlation  coefficient  becomes 


where: 


Xq  =  wavelength  corresponding  to  the  center  frequency 
Af  =  flat  bandwidth 

The  effect  of  bandwidth  on  the  correlation  coefficient  can  be  observed  in 
Figure  29  where  bandwidths  of  25%  and  50%  are  compared  with  the  single  frequency 
case  (zero  bandwidth). 

It  is  shown  in  Section  IV  that,  for  a  uniform  distribution  of  spatial  noise,  op¬ 
timization  of  the  signal -to -noise  ratio  is  synonymous  with  the  maximization  of  array 
directivity.  If  the  noise  power  is  uniformly  distributed  over  the  frequency  band  of 
interest,  this  corollary  holds  for  the  wideband  array.  Therefore,  the  effect  of  band¬ 
width  on  signal -to -noise  ratio  can  be  determined  by  applying  the  correlation  coefficients 
given  by  equation  66  to  the  conditions  of  maximum  directivity  presented  in  Section  IV. 
The  identity  equation  (Equation  29  in  Section  IV)  which  is  to  be  minimized  for  maximum 
directivity  is  expressed  as  follows  for  the  linear  wideband  array  of  isotropic  elements 
with  uniform  noise. 
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Figure  29.  The  Correlation  Coefficient  for  Band-Limited  White  Noise 
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Equation  67  can  be  minimized  by  solving  for  the  amplitude  coefficients  (amplitude 
weighting)  or  the  spacing  coefficients  (density  weighting)  as  in  Section  IV. 

The  effect  of  bandwidth  on  array  directivity  is  illustrated  in  Figures  30 
and  31  .  Figure  30  shews  the  optimized  directivity  as  a  function  of  amplitude 
control  for  bandwidths  of  0,  25%,  and  50%.  The  array  consists  of  five  elements  and 
the  equal  amplitude  cases  are  shown  as  dashed  lines  for  comparison.  The  five  ele¬ 
ment  array  is  optimized  as  a  function  of  element  spacing  in  Figure  31  with  the  equal 
spacing  case  added  for  comparison. 

The  optimized  directivity  response  curves  in  Figures  30  and  31  reveal 
several  significant  characteristics.  It  is  noted  that  the  maximum  directivity  is  reduced 
as  the  bandwidth  is  increased,  and  the  maximum  directivity  occurs  at  closer  element 
spacing  as  the  bandwidth  is  increased.  The  increase  in  array  directivity  available 
through  amplitude  control  (Figure  30  )  is  very  slight  regardless  of  element  spacing, 
while  the  optimized  element  spacing  (Figure  31  )  provides  a  considerable  increase 
for  large  element  spacing.  This  effect  is  the  result  of  spurious  lobe  suppression 
which  is  obtained  as  a  by-product  in  the  optimized  element  spacing  case. 

Since  amp1  ide  weighting  and  element  spacing  are  variables  in  the  optimiza¬ 
tion  of  array  directivity  in  the  wideband  case,  it  is  conceivable  to  utilize  both  variables 
simultaneously.  An  example  of  this  condition  is  presented  in  Table  5  .  The  conditions 
for  this  computation  are: 

a.  Five  element  array 

b.  Broad  side  fixed  beam 

c.  Isotropic  antenna  elements 

d.  No  mutual  coupling  between  elements 

e.  Uniform  noise  over  bandwidth  and  space 

1.  No  internal  uncorrelated  noise 

g.  Conjugate  matched  flat  receiver  response 

Corresponding  to  Figures  30  and  31 ,  it  is  noted  that  as  the  bandwidth  is  increased 
the  total  length  of  the  optimized  array  is  significantly  reduced.  While  the  element 
spacing  varies  considerably,  the  urnobtude  weighting  is  very  nearly  uniform.  From 
the  above  examples,  it  can  be  concluded  that 

a.  Increased  bandwidth  decreases  the  maximum  signal -to -noise  ratio 
in  an  array. 

b.  Optimum  performance  exceeds  the  uniform  amplitude  and  uniform 
spacing  ease  only  slightly,  except  at  large  element  spacing. 

c.  Optimized  performance  is  achieved  with  less  aperture  as  the  frequency 
band  increases. 

3.  SPACE  -  FREQUENCY  EQUIVALENCE 

In  the  wideband  array,  a  relationship  exists  between  element  spacing  and  band¬ 
width  which  can  be  used  as  a  design  parameter.  This  'oncept  was  originally  proposed 
by  Kock  and  Stone,  Reference  (21)  .and  later  adapted  to  wideband  arrays  by  Dausin,  et.al, 
Reference  (22).  The  application  of  space  frequency  equivalence  is  illustrated  in  Figure  32 . 
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Figure  30.  Maximized  Antenna  Directivity  by  Amplitude  Control 
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Figure  32.  Space-Frequency  Equivalence 
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The  location  of  any  array  element,  relative  to  the  center  line  of  the  array,  when  meas¬ 
ured  in  wavelengths,  is  a  function  of  frequency.  For  a  given  physical  location,  the  ele¬ 
ment  is  a  greater  number  of  wavelengths  removed  at  a  higher  frequency  (\  in  Figure  32) 

than  at  a  lower  frequency  (X^).  In  the  wideband  array,  this  variation  in  effective  dis¬ 
placement  is  referred  to  as  the  equivalent  location.  It  is  noted  that  as  the  physical  dis¬ 
placement  of  the  element  increases,  with  respect  to  the  array  centerline,  the  equivalent 
location  also  increases. 

The  utilization  of  this  concept  in  the  design  of  wideband  arrays  is  dependent  on 
the  ultimate  array  application.  To  illustrate  the  design  procedure,  consider  an  array 
with  an  overall  aperture  of  approximately  50Aq,  which  is  to  be  excited  with  25  spectral 

frequency  components  separated  by  0. 01  f  for  a  total  bandwidtn  of  25  percent. 

As  a  first  attempt  at  evaluating  the  space -frequency  equivalent  array,  the  ele¬ 
ments  are  arranged  so  that  the  equivalent  locations  just  touch  but  do  not  overlap.  For 
the  elements  near  the  array  center,  the  equivalent  locations  cannot  touch,  therefore 
these  elements  are  left  in  their  \q/2  locations.  The  resulting  array  contains  a  total  of 

thirty  elements.  The  computed  radiation  pattern  (at  t=  o)  is  shown  in  Figure  33  along 
with  the  pattern  of  the  same  array  operating  at  a  single  frequency  f  .  Since  complete 

evaluation  of  any  wideband  array  must  include  time  as  well  as  space  variations,  the  time 
response  for  several  values  of  u  is  shown  in  Figure  34  . 

Although  the  space-frequency  equivalence  has  the  effect  of  lowering  the  sidelobes 
in  the  far-out  region,  the  smearing  effect  shown  in  Figure  33  is  not  desirable  for  the 
purpose  of  producing  a  well-defined  major  lobe. 

It  is  noted  that  in  the  thirty  element  array,  the  first  six  elements  on  either  side 
of  the  center  line  remain  on  the  nominal  Xq/2  spacing.  In  Section  III,  a  space  taper  pro¬ 
cedure  is  described  wherein  the  sidelobes  near  the  major  lobe  can  be  reduced  through 
incremental  element  spacing.  Applying  this  principle  to  the  center  twelve  elements  of 
the  thirty  element  array  produces  the  patterns  shown  in  Figures  35  and  36  .  It  is 
noticed  that  a  relatively  uniform  sidelobe  level  is  maintained  over  the  entire  range  of 
angle  shown. 

The  purpose  of  the  above  example  is  to  illustrate  the  combination  of  array  de¬ 
sign  techniques  and  the  latitude  of  control  available  where  the  variables  are  not  con¬ 
fined.  Confinement  of  these  variables  is  dependent  on  application,  but  it  is  felt  that  the 
combination  of  design  techniques  can  be  effectively  used  in  many  practical  applications. 
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Frequency  Equivalence 


Figure  35.  Thirty  Element  Space  Tapered  Array 


SECTION  VI 


APPLICATIONS 


In  the  preceding  sections,  the  effects  of  array  geometry,  signal  allocation, 
frequency  bandwidth,  and  noise  distribution  on  the  response  of  a  linear  array  are  dis¬ 
cussed  in  detail.  This  investigation  has  revealed  several  array  characteristics 
which  could  be  applied  to  enhance  the  performance  of  particular  systems  utilizing 
array  antennas.  Several  design  and  functional  concepts  of  application  are  discussed  in 
this  section.  These  are  not  intended  to  be  complete,  but  rather  indicative  of  the  man¬ 
ner  in  which  the  observed  effects  may  be  used  to  complement  other  system  require¬ 
ments.  In  addition  to  the  applications  relating  to  linear  arrays,  examples  are  provided 
which  illustrate  the  effect,  on  signal -to -noise  ratio,  of  applying  the  principals  to  two 
dimensional  planar  arrays. 

1.  SIGNAL-TO-NOISE  RATIO  IN  DIRECTIVE  ARRAYS 

In  Section  IV  it  is  shown  that  the  signal-to-noise  ratio  is  directly  related  to 
array  directivity  in  a  uniform  noise  environment.  It  is  further  shown  (see  Figure  42) 
that  the  array  directivity  is  a  function  of  the  number  of  elements  employed  and  the 
spacing  between  them.  As  the  elements  become  closely  spaced  (less  than  X/2)  an  in¬ 
crease  in  directivity  (super -gain)  is  possible  through  critical  control  of  the  element 
amplitudes.  Although  this  does  not  represent  a  true  increase  in  array  gain,  it  can 
represent  an  increase  in  signal-to-noise  ratio  since  the  array  is  more  directive,  there¬ 
by  discriminating  against  its  noise  environment. 

If  a  particular  application  can  utilize  the  increased  directivity  in  the  array 
pattern,  and  the  demand  for  additional  elements  can  be  tolerated,  then  an  increase 
in  signal-to-noise  ratio  can  be  obtained  which  is  consistent  with  the  stability,  noise 
figure,  and  bandwidth  of  the  weighting  components.  For  example,  a  super  directive 
array  could  be  considered  in  a  radio  astronomy  application  where  the  array  is  required 
to  receive  only,  and  noise  and  directivity  are  of  paramount  importance.  Low  noise 
amplifiers  in  a  controlled  environment  would  be  required  to  establish  the  critical 
weighting  function. 

2.  SPECTRAL  OR  WAVEFORM  DETECTION 

It  is  shown  in  Section  IV  that  the  time  response  of  a  wideband  array  is  a  func¬ 
tion  of  the  angle  of  observation  with  respect  to  the  beam  maximum.  It  follows  there¬ 
fore,  that  the  observed  waveform  and  its  associated  spectral  content  are  also  dependent 
on  the  observation  angle.  The  only  angle  at  which  the  impressed  spectrum  or  wave¬ 
form  can  be  observed  is  boresight,  since  the  angular  position  of  the  beam  maximum 
is  not  a  function  of  frequency. 

The  spectral  response  which  is  observed  at  any  angle  other  than  boresight 
is  simply  the  antenna  response  to  each  of  the  spectral  components.  This  car  be  com¬ 
puted  for  any  regular  array  by  considering  the  spectral  components  of  the  impressed 
signal  separately  in  determining  the  array  response  as  a  function  of  the  angle.  The 
angular  spectrum  is  then  constructed  by  picking  the  response  at  a  particular  angle  for 
each  of  the  components  comprising  the  total  spectrum. 
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Since  the  impressed  waveform  and  its  associated  spectrum  are  accurately 
reproduced  only  in  the  major  lobe,  it  appears  that  using  the  impressed  waveform  as  a 
reference,  a  means  of  major  lobe  identification  can  be  established.  As  an  example, 
consider  the  pulsed  transmitter  excitation  of  a  uniform  array  depicted  in  Figure  37. 

The  boresight  response  is  shown  to  be  a  replica  of  the  impressed  signal  while  the 
sidelobe  response  is  distorted.  To  be  exact,  it  must  be  stated  that  the  input  is 
accurately  reproduced  only  in  the  center  of  the  major  lobe,  but  for  practical  purposes, 
the  shape  is  well  maintained  within  the  half-power  angle. 

In  a  one-way  propagation  path  the  transmitted  signal  would  be  collected  by  a 
receiving  antenna  and  processed  in  a  receiver  circuit  to  yield  the  desired  output.  This 
processing  might  take  the  form  illustrated  in  Figure  38  wherein  the  received  signal  is 
amplified  and  passed  through  a  filter  of  special  design.  If  the  filter  is  a  matched 
filter,  Reference  (23)  i.e.  a  conjugate  network,  maximum  output  can  only  be  obtained 
in  the  boresight  region.  The  filter  will  discriminate  against  the  distorted  waveform 
in  the  sidelobe  region  thereby  reducing  the  effect  of  sidelobe  radiation.  The  filter 
can  he  matched  to  the  total  waveform,  or  discrimination  can  be  achieved  on  a  spectral 
basis,  by  matching  to  the  specific  frequency  components. 

An  alternate  form  of  the  detection  technique  is  shown  in  Figure  39.  A 
sample  of  the  received  signal  is  processed  on  a  parallel  basis  and  used  to  control  a 
switch  which  completely  isolates  the  receiver  when  the  incoming  wave  does  not  match 
the  known  transmitter  function.  There  are  two  characteristics  of  parallel  signal 
processing  which  may  be  advantageous  to  some  systems.  When  the  switch  is  activated 
to  the  off  position,  complete  suppression  of  the  undesired  radiation  is  obtained  rather 
than  a  reduction  resulting  from  the  filtering  action.  The  parallel  processing  also  per¬ 
mits  the  received  signal  to  be  transferred  directly  to  the  receiver  without  modificatkm 
by  the  filter  which  may  be  desirable  in  some  applications. 

The  spectral  detection  technique,  as  herein  described,  presupposes  that,  (1) 
the  spectrum  of  the  transmitted  waveform  is  known,  and  (2)  that  it  does  not  vary.  If 
Lne  former  is  true  but  not  the  latter,  an  adaptive  filter  could  be  employed  which  is  ad¬ 
justed  to  variations  in  the  transmitter  spectrum. 

In  application,  the  two  way  path  of  a  radar  system  might  be  considered.  A 
replica  of  the  transmitted  spectrum  is  always  available.  If  the  radar  utilizes  a  pulsed 
transmitter,  the  spectrum  of  the  output  is  essentially  constant  being  determined  by  the 
pulse  characteristics  and  system  bandwidth.  If  the  target  is  assumed  to  approximate  a 
point  target,  the  received  signal  spectrum  will  be  essentially  that  transmitted  with  re¬ 
duced  amplitude.  The  spectral  detection  technique  could  then  be  employed  to  reduce 
the  spurious  response  of  sidelobe  and  gating  lobe  radiation 

In  essence,  the  processing  of  the  signal  spectrum  in  this  manner  is  a  self 
referencing  sidelobe  cancellation  technique.  The  assumption  which  would  limit  its 
usefulness  in  a  practical  radar  application  is  the  point  target  requirement.  Most  radar 
targets  are  not  point  targets  or  even  symmetrical  targets,  therefore,  the  observed 
spectrum  of  the  returned  signal  is  not  identical  to  that  transmitted.  A  probability 
criterion  could  be  established  for  specific  targets  but  the  effectiveness  of  the  detection 
scheme  would  be  reduced  with  the  possibility  of  lost  targets. 

3.  SECURE  COMMUNICATIONS 

An  adaptation  of  the  spectral  detection  technique  might  be  employed  to  secure 
a  communications  link  in  which  coded  information  is  to  be  relayed.  This  application, 
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Figure  37.  Transmitted  Waveform  and  Spectrum 
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Figure  38.  Direct  Signal  Processing 
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Figure  39.  Parallel  Signal  Processing 


80 


like  the  detection  concept,  is  based  on  the  fact  that  waveform  distortion  occurs  in  the 
far  field  of  the  array  for  all  angles  displaced  from  the  major  lobe. 

Assuming,  for  example,  that  the  information  to  be  transmitted  is  in  fhe  form 
of  short  pulses,  the  pulse  train  could  be  developed  to  produce  a  highly  distorted  re¬ 
sponse  in  the  sidelobe  region.  The  degree  of  distortion  is  dependent  on  the  total  band¬ 
width  required  for  transmission,  the  shape  of  the  pulses,  the  separation  between 
pulses,  and  the  characteristics  of  the  array. 

4.  ADAPTATION  IN  A  NON-UNIFORM  NOIP--  ENVIRONMENT 

One  of  the  conditions  assumed  in  the  study  of  signal-to-noise  optimization  is 
uniform  noise  distribution  over  all  space.  It  has  been  shown  that  with  this  confinement, 
optimum  signal-to-noise  ratio  occurs  for  near  uniform  weighting  and  equal  element 
spacing  of  approximately  300  degrees  in  the  narrow  band  case.  It  is  well  known  that 
uniform  noise  is  not  a  practical  noise  model,  even  in  the  "window"  which  occurs  be¬ 
tween  cosmic  noise  and  molecular  absorption  (1  GHz  to  10  GHz). 

Although  the  definition  of  a  practical  noise  model  is  beyond  the  scope  of  this 
project,  it  is  of  interest  to  determine  the  behavior  of  the  optimization  process  in  a 
specified  non-uniform  environment.  To  illustrate  a  possible  adaptive  control  tech¬ 
nique,  three  non-uniform  noise  models  are  considered  as  examples. 

a.  Uniform  noise  with  a  20  DB  increase  in  the  overall  sidelobe  region. 

b.  Uniform  noise  with  a  20  DB  increase  applied  only  in  the  far  out 
sidelobe  region. 

c.  Uniform  noise  with  a  10  DB  increase  in  the  major  lobe  and  near 
sidelobe  region. 

The  signal-to-noise  ratio  is  optimized  for  each  of  the  assumed  models  through  ampli¬ 
tude  control.  The  resulting  radiation  patterns  are  shown  in  Figures  40  through  42. 

It  is  noted  that  the  optimized  radiation  pattern  is  severely  modified  by  the  noise 
environment.  The  signal-to-noise  ratio  with  increased  sidelobe  noise  is  optimized  by 
maintaining  a  low  sidelobe  level  (Figure  40).  When  the  increase  in  noise  is  confined 
to  the  far  out  sideiobes  (Figure  41),  the  close  in  lobes  rise.  When  the  noise  is  in¬ 
creased  in  the  region  of  the  major  lobe  (Figure  42),  high  uniform  sideiobes  are  char¬ 
acteristic  of  the  optimized  pattern. 

This  variation  in  radiation  pattern,  as  a  given  array  is  optimized  for  maxi¬ 
mum  signal-to-noise  ratio  in  different  noise  environments,  suggests  the  possibility  of 
adapting  an  array  to  its  noise  environment  by  amplitude  weighting  of  the  antenna  ele¬ 
ments  . 


The  concept  of  noise  adaptation  could  be  implemented  in  a  number  of  ways. 
Consider,  for  example,  a  scanning  receiver  array  (mechanical,  or  electronic)  with  a 
variable  gain  amplifier  on  each  antenna  element.  With  the  g  in  of  each  amplifier  con¬ 
stant  (uniform  illumination)  the  array  could  be  scanned  to  determine  the  local  noise 
environment.  The  gain  of  the  element  amplifiers  could  then  be  adjusted  for  maximum 
signal-to-noise  ratio  as  described  above.  This  process  could  be  up-dated  at  any  time 
to  account  for  changes  in  the  environment  such  as,  antenna  orientation,  atmospheric 
changes,  or  the  presence  of  ECM  equipment. 
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Figure  40.  Optimizee  u 

Figure  40.  Optimized  S/N  for  a  20  DB  Noise  Increase  in  the  Sidelobe  Region 
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Figure  41.  Optimized  S./N  for  a  20  DB  Noise  Increase  in  the  Far-out  Sidelobe  Region 
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Figure  42.  Optimized  S/N  for  a  10  FB  Noise  Increase  Near  Boresight 


5.  APPLICATION  TO  PLANAR  ARRAYS 


The  primary  area  of  Investigation  in  this  study  is  the  linear  array.  A  two 
dimensional  extension  of  the  generalized  cross-correlation  coefficient  (Appendix  C) 
indicates  that  the  conditions  for  optimizing  the  planar  are  of  the  same  form  as  those 
encountered  in  the  analysis  of  the  linear  array.  The  complicating  factor  in  the  planar 
array  analysis,  is  the  two  dimensional  element  spacing  factor  which  permits  the  array 
to  assume  an  unconfined  geometrical  form.  For  comparison,  two  common  element 
arrangements  are  considered:  (1)  square  spacing  and  (2)  triangular  spacing.  The  op¬ 
timized  directivity  (S/N  in  the  uniform  noise  environment)  is  snown  in  Figure  43  for 
the  narrow  frequency  band  and  in  Figure  44  for  a  50  percent  bandwidth.  In  each  case, 
the  planar  array  is  comprised  of  25  equal  amplitude  elements. 

The  results  of  the  planar  array  closely  parallel  those  obtained  for  the  linear 
array  (See  Figures  30  and  31).  It  is  noted  that  while  maximum  directivity  is  obtained 
with  square  element  spacing,  the  triangular  spacing  produces  higher  directivity  for 
closely  spaced  elements.  This  brief  investigation  illustrates  that  the  techniques 
developed  for  the  optimization  of  signai-to-noise  ratio  in  linear  arrays  apply  also  to 
the  planar  array  although  the  computational  procedure  is  more  complex. 
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gure  43.  Optimized  Planar  Array  -  Narrow  Band 


Figure  44.  Optimized  Planar  Array  -  Wideband 


SECTION  VII 


COMPUTATIONAL  TECHNIQUES 


During  the  course  of  investigation  covered  by  this  report,  several  significant 
mathematical  techniques  evolved,  which  are  described  in  this  section.  Each  technique 
is  described  qualitatively:  the  details  of  analysis  of  the  various  methods  appear  in  the 
appendix  Examples  are  included  where  appropriate. 

1.  A  METHOD  FOR  CONTROLLING  PATTERN  SIDELOBE  LEVEL  THROUGH 
SPACE  AND  AMPLITUDE  TAPERING 

This  method  evolved  from  a  method  by  Harrington,  Reference  (24).  Basically, 
the  method  is  an  iterative  technique  using  the  finite  difference  approach.  It  was 
originally  intended  to  be  applied  to  space  tapered  arrays,  but  during  the  development 
of  the  method,  it  became  clear  that  the  technique  is  general  and  therefore  is  applicable 
to  amplitude  tapered  arrays. 

The  normalized  far  field  pattern,  for  an  array  of  2N  elements  equally  spaced  is 
given  by, 


E<“>  3  lE  Z  Ar.  cos  ^(n-0.5)7ruj 


where 


^  (sin  9  -  sin  9q) 


d  =  is  the  element  spacing 

X  =  is  the  wavelength  at  the  operating  frequency 

9  =  is  the  observation  angle  from  boresight 

9  =  is  the  array  steering  angle 

th 

An  =  is  the  amplitude  of  the  n  element 

Applying  finite  differences  to  both  sides  of  equation  .68 
N 

AE  AAn  cos  [^(n-0.5)7ruj 

y  a  n=l 

X  ,  n 


The  far  field  pattern  is  then  described  by, 
Erf(u)  =  E  (u)  +  AE 
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Implementation  of  this  method  consists  of  introducing  small  changes  in  the  ampli¬ 
tude  coefficients  of  the  normalized  far-field  pattern  equation  which  are  proportional 
to  the  sum  of  the  amplitudes  of  an  arbitrary  number  of  impulse  functions,  each  of 
which  is  multiplied  by  a  cosine  function  whose  argun  ent  involves  the  points  at  which 
the  impulse  functions  are  applied. 

The  ,-mpulse  functions  are  applied  at  the  sidelobe  maximum  of  the  power  pattern 
and  the  number  of  such  functions  used  depends  upon  how  many  sidelobes  require  re¬ 
duction.  The  degree  of  reduction  is  determined  by  the  absolute  difference  in  the  de¬ 
sired  and  actual  sidelobe  peaks.  Reasonable  sidelobe  control  can  be  obtained  by  re¬ 
ducing  only  one  sidelobe  at  a  time,  although  some  reduction  in  computation  time  can 
be  achieved  by  applying  impulses  at  all  of  the  sidelobe  maxima  for  the  range 
0°  <  9  <  90°.  These  amplitude  deviations  are  then  added  algebraically  to  the  original 
amplitude  coefficients,  which  can  be  arbitrary.  The  resulting  voltage  pattern  is  com¬ 
puted  using  these  new  coefficients  and  the  cycle  is  repeated.  This  iteration  is  con¬ 
tinued  until  the  resultant  pattern  approximates  the  desired  pattern  within  a  prescribed 
error 

Application  of  the  iterative  technique  is  simplified  through  the  application  of  a 
side  lobe  deviation  function 

K 

AAn  =  ~F  2^  ak  C0S  5)  ^k1  (71) 

k=l 

where 

ak =  E<uk)  -  Ed  ‘V  m 

and  F  is  a  convergence  factor,  u^  is  the  point  at  which  the  impulse  function  is  applied, 
and  K  is  the  number  of  sidelobesK  to  be  reduced. 

Equations  71  and  72  can  be  directly  applied  to  a  computer  program  for  solution. 

The  programming  procedure  desc~  'bed  in  the  algorithm  of  Table  6. 

To  ensure  convergence,  it  may  be  necessary  to  modify  the  amplitude  deviations 
by  a  scaling  factor.  This  scaling  factor  can  be  determined  by  trial  and  error  until 
convergence  of  the  method  is  obtained. 

The  derivation  of  equations  71  and  72  appears  in  Appendix  A. 

2.  AN  EXACT  METHOD  OF  DETERMINING  THE  ARRAY  OUTPUT  WAVEFORM 

FOR  ANY  SPECIFIED  INPUT  WAVEFORM 

This  technique  is  based  on  the  application  of  the  array  impulse  response  to  the 
far-field  pattern  equation  containing  the  frequency  characteristics  of  the  antenna  for  a 
given  observation  angle.  The  inverse  Fourier  Transform  is  then  applied  to  the  product 
of  the  array  input  spectral  density  function  and  the  far  field  pattern  function.  The  re¬ 
sult  is  the  impulse  response  of  the  array  if  the  array  input  spectral  density  function  is 
considered  to  be  a  constant.  Simplification  of  the  impulse  response  can  be  effected  by 
sampling  the  function  at  specific  times.  The  convolution  integral  is  then  applied  to 
the  array  impulse  response  and  the  specified  input  time  waveform,  determining  the 
output  time  waveform.  The  waveform  thus  obtained  is  described  by 
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TABLE  6  .  ALGORITHM  FOR  ITERATIVE  PATTERN  SYNTHESIS 


1.  Initiate  the  iterative  program  vith  arbitrary  array  coefficients. 

2.  Compute  the  normalized  voltage  pattern. 

3.  Tabulate  the  sidelobe  peaks,  location  of  the  peaks,  and  number  of  sidelobes 
within  the  desired  control  region. 

4.  Compute  the  impulse  magnitudes,  a^,  at  the  sidelobe  peaks,  u^. 

5.  Compute  the  amplitude  (or  spacing)  correction  for  each  element  using 
equations  71  and  72. 

6.  Apply  the  correction  to  each  of  the  amplitudes  (or  spacing). 

7.  Go  to  step  2  if  last  pattern  computed  is  not  sufficiently  approximated  by  the 
desired  pattern,  otherwise  proceed  to  the  next  step. 

8.  List  the  element  amplitudes  (or  spacing). 

9.  Compute  the  resulting  pattern. 
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-itrnu , 

V (u,  t)  =]TAn  e  v  T(t-tn) 
n 


(73) 


where 


4n  =  rr  <Vu> 

o 

is  the  amplitude  of  the  nth  element 
T(t-tn’i  is  the  specified  input  time  waveform 
t  is  the  nth  samoling  time 
u^  is  the  phase  steering  angle 
ut  is  the  real-time  steering  angle 

The  above  equation  is  a  simp'a  finite  summation  which  is  easily  applied.  This 
method  is  applicable  to  any  discrete  array. 

As  an  example  of  the  application  of  ''his  method,  consider  a  simple  CW  signal 
expressed  as, 

trf  t 

(74) 


j2irf  t 


T(t)  =  e 

From  equation  73  the  waveform  is  described  by 


j2zfi 


V(u,  t)  -  e  °  ]TAne 


^(u-UfV 


(75) 


Equation  75  is  the  general  time,  space  antenna  pattern  of  an  array,  subject  to 
both  phase  and  real-time  steering.  The  time  variation, 
j2rfol 
e 

represents  the  carrier  frequency. 

As  another  example,  consider  a  single  linear  FM  pulse  defined  by 


t  •>  ^ 

h)e 

The  array  response  to  this  signal  is 

) 

V  A 


T(t)  =  RECT.  (  --  )  e 
1 


9  / ,  nu 

j(2*ft+  ;;  t  )  (  ’  2f 

V(u,  t)  =  e  J  v"  A_  RECT.  1  ° 


T  T 

1  1 
<  t  <  — i- 
2  -  -  2 


inuV 


(76) 


(77) 


92 


where 


V 


+ 


Mnu 

2 

8 1 


o 


and  u  =  u  =  0. 

«w  t 


3.  MAXIMIZATION  TECHNIQUES 


Several  methods  of  mathematical  processing  are  available  to  the  problem  of 
maximizing  array  directivity  or  signal-to-noise  ratio  with  respect  to  element  ampli¬ 
tude  control  or  element  spacing.  All  of  these  methods  can  be  generally  catagorized 
as  either  an  iteration  process  or  a  matrix  inversion  technique.  The  choice  of  method 
is  entirely  dependent  on  the  parameters  of  the  problem.  The  iterative  technique  is 
often  more  accurate  and  can  be  applied  to  large  arrays,  but  for  some  conditions  fails 
to  converge.  The  matrix  inversion  method,  on  the  other  hand,  is  difficult  to  imple¬ 
ment  for  large  arrays  and  is  subject  to  round  off  errors.  A  combination  of  these 
techniques  can  be  used  to  advantage  in  some  specialized  cases,  but  in  general  the 
combination  incurs  some  of  tne  disadvantages  of  each  method. 

a.  Iterative  Processes 

The  iterative  technique  of  maximizing  a  function  follows  the  process  described 
in  Table  7  .  In  beneral  this  process  is  applied  to  the  first  or  second 
derivitive  of  the  function  to  be  maximized.  An  example  of  this  process  is 
described  in  Section  IV,  4  for  determining  the  maximum  array  directivity 
from  control  of  the  applied  amplitude  coefficients. 

The  primary  advantage  of  the  iterative  process  is  the  simplicity  of  form. 

A  large  number  of  terms  can  be  handled  in  computation  with  minimum 
capacity.  The  computer  time  required  to  reach  a  solution  may  be  considerable 
if  the  solution  does  not  converge  rapidly,  but  for  a  function  with  many  terms 
this  may  be  the  only  form  of  solution.  The  major  disadvantage  of  the.  iterative 
technique  is  the  fact  that  optimum  performance  is  not  assured.  The  calculation 
may  converge  about  any  maximum  value,  not  necessarily  the  optimum  or.e. 
Parallel  computation  can  be  used  as  a  check  on  optimization. 

b.  Matrix  Inversion 

The  problem  of  matrix  inversion  has  been  investigated  for  many  years  and 
several  excellent  forms  of  solution  have  been  obtained,  Reference  (25).  Since 
the  problems  encountered  in  this  study  are  specific  rather  than  general,  only 
a  form  utilizing  Gaussian  elimination  has  been  applied. 

The  matrix  inversion  technique  described  has  proven  effective  in  obtaining 
the  inverse  of  a  non-singular  matrix.  The  only  prerequisite  for  the  matrix 
is  that  none  of  its  diagonal  elements  be  zero. 
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TABLE  7  .  MAXIMIZATION  ALGORITHM 

1.  Assume  an arbj.trar' set  of  amplitude  coefficients;  A^,  A^,  . . . ,  A^. 

2.  Solve  the  iteration  equation  A^  =  QA^  +  R  for  all  values  of  k: 

k  =  1,  2,  . . . . ,  N.  In  general,  the  factors  Q  and  R  are  functions  of  n  and  An. 
(See  Section  IV. 4  of  this  report.) 

3.  If  the  new  coefficients  agree  within  a  prescribed  error  with  the  old  coefficients 
terminate  the  iteration.  Otherwise,  go  to  the  next  step. 

4.  Replace  the  old  set  of  amplitude  coefficients  with  the  new  set. 

5.  Go  to  step  2. 
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Let  the  nxn  matrix  [A]  be  denoted 


The  basic  matrix  is  then  augmented  with  an  nxn  identity  matrix  [I]  to  obtain 
a  modified  matrix,  [A/I] 


By  a  sequence  of  elementary  row  operations  upon  A/I,  an  attempt  is  made  to 
transform  the  left-half  of  [A/I]  into  an  identity  matrix.  If  this  is  possible 
then  it  can  be  shown  that  the  right  half  of  79  is  then  [A ]“*. 


The  first  portion  of  the  procedure  is  to  obtain  zeros  in  the  off-diagonal  terms 
of  A  and  unity  along  the  diagonal.  This  procedure  is  best  illustrated  by  means 
of  an  arbitrary  example.  Consider  the  following  3x3  matrix. 

'  1  -2  3 

[Ai  =-3  7  12 

4  -4  -2 

Augmenting  this  matrix  with  the  identity  matrix  [I] 


(81) 
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It  can  be  shown  that  a  valid  operation  on  81  is  to  acid  to  any  row  (column) 
any  multiple  of  another  row  (column).  Thus  in  81  ,  to  obtain  a  zero  in  row 
2  and  column  1,  we  can  add  row  1  multiplied  by  3  to  obtain: 


1 

-2 

3  : 

L  0 

o' 

[A/I]'  = 

0 

1 

21  : 

1  1 

0 

(82) 

L4 

-4 

-2  1 

[)  0 

1 

The  next  step  is 
which  yields: 

to  add 

the  first  row  multiplied  by  a  negative  4  to  row  3; 

l 

-2 

3 

1  0 

o' 

[A/I]"  = 

0 

1 

21 

3  1 

0 

(83) 

0 

4 

-14 

-4  0 

1 

The  elimination  of  the  element  in  the  third  row  and  second  column  is  achieved 
by  adding  row  2  multiplied  by  -4  to  row  3; 


[A/I] 


-2310 
1  21  3  1 

0  -98  -16  -4 


0 

0 

1 


(84) 


Repeating  this  procedure  on 

the  upper 

triangular 

part  of  [A]  yields, 

"l 

0 

0 

-0.349 

0.163 

0.459  ' 

[1/A‘V 

0 

1 

0 

-0.43 

0. 143 

0.  214 

_0 

0 

1 

0.163 

0. 0408 

-0.0102 

(85) 


which  is  the  desired  inverted  matrix.  A  problem  consisting  of  a  set  of 
simultaneous  equations  can  be  solved  by  using  the  above  method  or  a  partial 
Gaussian  elimination  procedure  may  be  incorporated. 

— •  — * 

The  partial  procedure  may  be  applied  to  an  nxn  matrix  (real)  and  X  and  Y 
which  are  two  nxl  column  matrices 

[A]  X  =  (86) 


The  procedure  is  to  obtain  the  form  of  equation  86  above.  If  the  left-half  of 
the  matrix  in  84  is  identified  by  A*  and  the  right  half  as  I*  then  equation 
becomes 

[A]*  [If  (87) 


In  the  example  cited  previously,  this  expression  becomes, 
^  =  -  ■k  <-16Yr4Y2+Y3> 

X2'3VY2-21X3 


X. 


=  Y1+2X2-3X3 


(88) 
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Thus,  the  solution  is  apparent  since  Y  is  specified  and  X  *  =  [X 

The  solution  vector  of  a  set  of  linear  equations  is  obtained  faster  on  a  digital 
computer  by  this  method  of  partial  elimination  than  by  calculating  the  inverse 
using  full  Gaussian  elimination. 

Combination 


A  combination  of  techniques  may  be  applicable  in  the  solution  of  specific 
problems.  Since  there  is  an  infinite  number  of  combinations  possible,  the 
technique  is  described  as  an  example  wherein  the  gain  of  an  antenna  array  is 
to  be  optimized.  A  form  of  the  gain  function  can  be  expressed  as  the  ratio 
of  two  bilinear  forms, 


G(A) 


and 


J3 


I  z 

n  m 

A  ■*  a 
m  n,  m 

T  T  A  *  p 
^  m  n.m 

n  m 

27rd  , 

sin 

—  (n-m) 

n,  m 


ffd  , 

T  (n'm> 


(89) 


(90) 


“i}  =  1  t.j -i.a.3 . N 

Obviously,  [»*j  is  a  singular  matrix  but  this  does  not  present  a  problem  in  the 
development. 


Applying  a  theorem  of  Gantmacher,  Reference  (26),  the  condition  for  which 
G(A)  is  maximum  is  given  by, 

<91) 

A  is  a  nxl  column  vector 


A 


(92) 


Substituting  equation  92  inr.o  the  vector  form  of  equation  89 

G(A)  -  \ 


max 


(93) 
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The  characteristic  equation  corresponding  to  91  can  be  shown  to  be 

I  [“1_Xpiaxl/3]l  =°  <94> 

where  the  bars  denote  the  deterrcinent  of  the  enclosed  expression. 
Premultiplying  both  sides  of  91  by  the  inverse  of  [j3]  and  transposing  terms 

{ipfVi-  w1]}^0  (95) 

Defining  a  new  matrix  [y], 

{M-XmaxlI]}  A=°  (96) 

where  y..  =  ^  P  •  K  is  important  that  [of]  be  premultiplied  on  the  left  by 

j 

[/3]  *  and  not  on  the  right.  Thus  all  columns  of  [y]  are  identical.  This 
matrix  [y]  is  valid  only  for  the  single  frequency  case,  but  with  a  modifica¬ 
tion,  a  similar  matrix  can  be  constructed  that  accounts  for  the  wideband 
case.  In  either  case,  the  computational  procedure  is  the  same. 

In  computation  the  inverse  of  [/?]  is  premultiplied  by  [of]  to  obtain  [yj .  The 
coefficients  of  [y]  are  then  obtained  by  an  iterative  process. 

-W0) 

Starting  with  an  initial  vector  X'  '  the  components  of  which  are  arbitrary,  the 
iteration  is  performed  according  to  the  following  algorithm. 

(1)  Let  n  =  0 

(2)  Compute  a  new  vector  X  ^n+1^  =  [y]  X'1^ 

(3)  Determine  the  value  of  the  element  of  X  that  has  the  largest 
absolute  value;  call  this  value  X  . 

(4)  Increment  n  by  one. 

(5)  If  X  -A  ,  <  e  or  X  ^  -  X  ^  <  5  exit  from  algorithm,  otherwise 

n  n-i  —  — 

go  to  step  2  above. 

Normally  only  two  or  three  iterations  are  needed  before  convergence  is 
obtained. 
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APPENDIX  A 

DERIVATION  OF  THE  PATTERN  SYNTHESIS  TECHNIQUE 

This  appendix  derives  the  appropriate  equations  and  describes  the  iterative  syn¬ 
thesis  technique  for  an  amplitude  tapered  array.  The  procedure  for  space-tapered  arrays 
is  similar  but  more  involved. 

Equations  68  ,  ,  69  ,  and  70  in  Section  VII  form  the  starting  point  for  the 
development  of  the  method.  Equation  70  can  also  be  written  in  another  form. 


Ed  (u)  =  — | -  £  (\+  AAr)cos  [  (n-.  5)  tt  u] 

E  A  n=1 


(A  -1) 


The  problem  is  to  determine  AE  so  that  the  desired  pattern  can  be  achieved. 
Some  consideration  must  be  given  as  to  how  the  desired  pattern  should  be  specified.  The 
antenna  properties  of  beamwidth  and  illumination  efficiency  are  primarily  a  matter  of 
a  tradeoff  with  the  desired  sidelobe  level  and  number  of  elements,  but  are  not  important 
here.  The  important  feature  of  an  antenna  pattern  is  the  envelope  of  the  sidelobe  peaks. 
It  is  natural,  then,  to  define  the  desired  pattern,  E^(u),  in  terms  of  the  sidelobe  level 
envelope.  The  phase  of  a  given  sidelobe  is  not  important  insofar  as  the  sidelobe  peak 
envelope  is  concerned  but  must  be  considered  when  used  with  equation  (A-l).  This 
amounts  to  assigning  a  positive  or  negative  sign  to  E^u),  depending  upon  the  sign  of 
E(u)  +  AE,  and  presents  no  practical  difficulty. 

Using  equations  69  and  70  from  Section  VIE  and  rearranging  terms: 


N  M 

c  /  N  \ 

Lt  A An  cos  [  (n-.  5}  tt  u]  =  (  £J  A  )  (E.(u)  -  E(u) 
n=l  \  n=l  / 


(A -2) 


The  left  hand  term  of  (A-2)  is  recognized  as  a  finite  Fourier  series  which  has  a  period  of 
-1  <  u  <  1.  Since  it  is  a  cosine  series  which  indicates  that  it  is  an  even  function,  and  the 
right  hand  term  is  also  an  even  function,  the  coefficients,  A  A  ,  can  be  found  from  the 
half-range  expansion 


[  E^(u)  -  E(u))  cos  [  (n-.  5)  tt  u)  du 


(A-3) 


n=l ,  2, 


The  value  of  this  integral  is  highly  dependant  upon  the  assumed  behavior  of  the 
desired  pattern,  E^(u).  Since  the  only  feature  of  interest  in  the  patterns  is  the  peak  side¬ 
lobe  level,  the  exact  behavior  of  E<j(u)  can  be  quite  arbitrary.  The  integral  can  be  evalu¬ 
ated  at  the  points  of  concern—  namely,  the  sidelobe  peaks  —  by  using  the  Dirac  delta 
function.  Thus, 


/  N  x  1 

2^  L  A nJ  J  !  Ed(u)  -  E(u)  ]  6  (u-uR)  cos 


[  (n-.  5)  tt  ]  du 


(A -4) 
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where  u,  are  the  locations  of  the  siuelobe  peaks  which  are  to  be  altered.  Performing  the 
indicates  integration  of  (A -4)  yields 


AA 

n 

/  N 
-  2  E 

'  A  )  E  a,  cos  [  (n-.  5)  tt  u,  ] 

(A -5) 

v  n=l 

k=l  K 

where 

ak 

-  E<\> 

-  Ed 

(A -6) 

Equation  (A -5)  may  be  written  in  the  form 

K 

AAn  =  F  Yj  ak  cos  [(n-.  5)tt  i^]  (A -7) 

k=l 


which  is  identical  to  equation  71  in  Section  VII. 

The  only  restraint  on  equation  (A-7)  is  concerned  with  the  F  factor.  For  a 
specified  voltage  pattern,  it  is  only  required  that  F  be  small  enough  so  that  the  iteration 
process  converges. 
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APPENDIX  B 
NOISE  ANALYSIS 


The  noise  power  which  appears  at  the  terminals  of  a  circuit  is  a  function  of 
the  circuit  parameters  and  environment.  Noise,  and  its  effect  on  the  performance  of 
a  system,  is  specified  in  a  variety  of  terms.  The  purpose  of  this  appendix  is  to  define 
the  terms  and  to  provide  the  relationships  which  can  be  used  in  circuit  analysis.  Ac¬ 
tive  and  passive  circuits  are  included  with  examples  to  illustrate  the  significance  of 
component  placement. 

1,  DEFINITIONS  AND  NOMENCLATURE 

The  following  terms  are  defined  and  assigned  symbols  for  future  reference. 

a.  Noise  Power  (N)  -  The  maximum  power  which  can  be  transferred 
from  a  resistor  of  value  R,  operating  at  temperature  T,  to  a  re¬ 
sistive  load  of  value  R. 


N  =  kTB 

where  k  =  Boltzmann's  Constant 

.-25 


(3-1) 


(1.38  x  10  joule/degree  K) 

T  =  Noise  Temperature  -  °K 
B  =  Bandwidth 

b.  Bandwidth  (B)  -  The  frequency  band  over  which  noise  is  considered. 


00 


B  = 


|H(f0)l2 
where: 


/  |H(f)|“  df 


(B-2) 


c. 


H(f>  =  frequency  response 
H(f  J  =  maximum  response 

(usually  center  of  band) 

Signal-to-Noise  Ration  (S/N)  -  The  usable  power  of  a  signal  divided 
by  the  noise  power  (N) 

S/N  .  jly  (B-3) 

d.  Noise  Figure  (F)  -  The  ratio  of  the  signal-to-noise  ratio  at  the 
input  of  a  circuit  to  that  appearing  at  the  output  of  the  circuit. 

(S/N)  in 

F  *  ,B-4) 


2.  DERIVATION 

Using  the  noise  figure  (F)  as  a  basic  quality  factor,  and  noise  temperature 
(T)  as  a  basic  circuit  parameter,  the  following  relationships  can  be  established. 
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Consider  a  basic  network  with  gain  G  and  operating  at  the  input  noise  temper¬ 
ature  T. 


kTB 


(S/N)  =  GS. 

'  '  'out  in 


kTBG  +  AN 


AN  =  noise  power  introduced 
by  the  ci  xuit 


From  eauation  (B-4) 

„  _  kTBGH-AN  _  AN 

kTBG  kTBG 


(B-5) 


The  noise  temperature  (T)  has  been  standardized  to  an  ambient  temperature 
of  290°K  (approximately  63°F).  This  specific  noise  temperature  is  designated  by  the 
symbol  T,.  If  the  above  network  is  operating  in  a  standard  environment,  equation 
(B-5)  becomes 


F 


1  + 


AN 

kT0BG 


(B-6) 


AN,  the  internal  noise  of  the  circuit,  can  now  be  defined  in  terms  of  the  effective 
noise  temperature  of  the  circuit  referred  to  the  input  terminals  (Te). 

AN  =  kTeBG  (B-7) 


Substituting  (B-7)  into  (B-6)  the  circuit  noise  figure  becomes 

T 

F  =  1  +  Tjr  (B-8) 

T 

0 

The  term  is  sometimes  referred  to  as  excess  noise  temperature  since  it  repre- 

o 

sents  a  deviation  from  the  standard  ambient  temperature  (T0). 


3.  CASCADED  CIRCUITS 

In  any  electronic  system,  the  overall  noise  figure  must  be  determined  from 
the  components  comprising  the  system.  Since  these  components  are  usually  cascaded, 
the  noise  interaction  must  be  considered. 

With  two  circuits  operating  in  tandem,  the  input  signai-to-noise  ratio  follows 
from  (B-3). 


"IRE  Dictionary  of  Electronics  Terms  and  Symbols"  Institute  of  Radio  Engineers, 
New  York  1961. 
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The  signal -to-noise  ratio  out  of  the  first  circuit  and  into  the  second  follows  from  the 
single  circuit  case. 


N2  “  kT0BG1  +  AN 

where:  AN^  =•  k'T^BG^ 

therefore:  S2  G1  Si 

=  kBG1  (T„  +  T“) 

and  the  signal-to-noise  ratio  at  the  output  of  the  second  circuit  becomes 

'  G1G2^1  Q181 
(S/N)out  =  kBG]G2{T0  +  Tx)  +  AN"  =  kBG  ^  (T  0  +  Tj)  +  kTg  B  (B-10) 

since  ANg  =  kTgBGg 


The  overall  noise  figure  (F0)  is  found  by  taking  the  ratio  of  equations  (B-9) 
and  (B-10), 

S. 


1 


F„  = 


kT0B 


(t0  +  t1)G1  +  t2 


G1S1 


k(T0  +  T^BGj  +  kTgB 


T„G, 


or 


F°  -  (!  +  T  )  +  t  G 


but  from  (B-S) 


F  =1+7— 

1  T0 


°  1 


and  F2  =  i  +  Tf- 


therefore: 


F0  =  Fj  + 


<F2 


(B-ll) 


If  this  basic  two  circuit  analysis  is  extended  to  the  general  n-circuit  case, 
the  overall  noise  figure  can  be  determined. 


(F9  ~1)  (F„  -1) 


F  =  F  + _ — 

0  1  G, 


G1G2 


+ - + 


(F  -1) 
'  n  ' 


Gl'  °2  Gn-1 


(B-12) 


or  by  substituting  from  equation  (B-8) 

Te  =  (F 0  -1)  T0 

Tg  (overall)  =  Tx  +  Tg  +  Tg  + — +  Tr 


(B- 13) 


G1G2 


G1’G2’  Gn-1 
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4. 


LOSSY  CIRCUITS 


In  the  foregoing  analysis,  it  was  assumed  that  the  circuits  are  active  with  a 
value  of  gain  (G)  greater  than  unity.  A  passive  lossy  component,  such  as  a  transmis¬ 
sion  line  section,  has  a  gain  of  less  than  unity  or  a  loss  factor  (L)  greater  than  unity 
(L  =  i/G).  As  in  the  active  component  case,  the  lossy  component  has  an  effective 
noise  temperature  (T  ).  Assuming  that  the  ambient  temperature  of  the  lossy  network 
is  T,  the  effective  noise  temperature  of  the  circuit  is  related  to  T  and  L.  The  total 
noise  out  of  the  circuit  is 

N  =  kT  B  =  +  AN  (B- 14) 

out  e  L  v  ’ 

With  no  internal  active  components,  the  noise  in  the  passive  circuit  comes  only  from 
the  loss,  therefore: 

AN  -  kT0B  -  kT^B  =  kT0B(l-  ^  )  (B-15) 

Substituting  (B-15)  into  the  basic  noise  figure  equation  (B-6),  the  overall  noise  figure 
of  the  lossy  component  becomes 

Fe  =  1  +  =  1+  L(l-T>  =  L  (B-16) 

Lossy  components  can  be  considered  in  cascade  with  active  components  by 
applying  equation  (B-16)  to  the  generalized  equation  (B-12)  for  the  appropriate  net¬ 
works. 

5,  COMBINATION  NETWORK  -  ACTIVE  AND  PASSIVE  COMPONENTS 

Using  the  cascaded  circuit  approach  outlined  above,  an  active  circuit  with 
gain  (G)  and  a  passive  circuit  with  loss  (L)  can  be  considered  as  two  discrete  cases. 

Case  I  -  The  lossy  component  (F  )  is  placed  in  front  of  the  active  circuit 
(F  )  in  the  signal  flow  path.  From  equation  (B-ll)  we  have 

a 

F0  =  F  +  (F  -1)  L  since  G  =  l/L 

0  SL 

and  from  (B-16)  Fa  =  L 

therefore  F0  =  L  +  (F  -1)  L  =  LF  (B-17) 

a  a 

Case  II  -  The  active  component  precedes  the  lossy  component  in  the  signal 
flow  path.  Again  from  equation  (B-ll) 

F0  =  F  +  t  -1  and  F  =  L  (B-18) 

a  0 _  0 

G 

The  general  noise  figures  for  the  two  cases  are  given  by  equations  (B-17)  and  (E-1S). 
To  realize  the  significance  of  component  location,  three  conditions  of  the  gain-loss 
relationship  are  considered. 
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Condition  A  -  G  »  L 

Case  I  -  F0  =  LF 

SL 

Case  IT  -  F0  =  F 
0  a 

For  any  value  of  loss  other  than  the  trivial  lossless  case,  it  is  advantageous  to  place 
the  amplification  ahead  of  the  loss. 

Condition  B  -  L  »  G 

Case  I  -  F0  =  LF 

a 

Case  II  -  Fe  =  F  +L 
e  a 

In  this  improbable  condition,  it  is  still  desirable  to  place  the  amplification  first  since 
the  product  of  two  numbers  greater  than  unity  is  generally  greater  than  or  equal  to 
their  sum. 

Condition  C  -  L=G 

Case  I  -  F„  =  LF 

a 

Case  II  -  F0  =  F  + 

a  Li 

The  loss  term  can  vary  in  magnitude  from  zero  to  unity  (as  a  function  of  L)  therefore 
Case  II  is  still  preferable  in  the  location  of  components. 
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A  PPENDIX  C 

CORRELATION  COEFFICIENTS 


where 


The  correlation  coefficient  for  a  linear  discrete  array  is  defined  as: 

1  .2tt  7r  J  — r ~  cos  0 

tt  J  J  /  I  (0,  P)  G  (6,  p)  W  (f)  e  sin6d6dpdf 

=  47r  -00  0  0 


/  W(f)df 


(C-l) 


G (0,0)  -  the  element  power  pattern,  (normalized  to  unity) 

1(6,  p)  =  the  noise  power  distribution  over  space,  (normalized  to  unity) 
W(f)  =  the  power  spectral  density 

The  S/N  ratio  of  the  array  output  for  a  bores  ight  beam  is 


N  N 

y.  V]  A  A  R 
*-*  Z-1  n  m  n;  m 

n=l  m=l _ 

/  N  %  2 


/  N 

E,  A" 
\n-l 


(C-2) 


where  the  peak  signal  is  taken  as  unity.  A  consequence  of  normalizing  the  peak  sig¬ 
nal  to  unity  is  that  the  restriction  R-l  at  d=0  must  be  applied. 

If  the  noise  power  distribution  is  uniform,  i.  e. ,  1(6,  p)  =  1  and  only  a  single 
frequency  is  considered,  theS/N  of  the  array  is  identical  to  the  directivity  of  the 
array. 

The  problem  of  determ inmg  the  correlation  coefficient  for  different  models 
is  primarily  that  of  carrying  out  the  integration  indicated  in  equation  (C-l). 

For  a  single  or  multifrequency  spectrum  we  have 


W(f)  =  H  Ar  5  (f-fk) 
k 


(C-3) 


r  W(f)  df  =  C  A 
_00  V  K 

for  any  element  and  noise  powei  distribution. 


(C-4) 


For  the  case  of  an  isot  pic  element,  single  frequency,  and  isotropic  noise, 
equation  (C-l)  reduces  to: 


•  271-  f  d 

i  2ir  7r  J  - — 

R(d)  =  t~  J  /  e  c  sin  0  d  3  d  0 

0  o 

Using  a  change  of  variable 


(C-5) 


7T  f  d 

u  =  - - —  cos  9 


du  = 


7T  f  d 
O 


sin  9 


(C-6) 


the  correlation  coefficient  assumes  the  form: 


R(d)  = 


(  2tt  d  \ 

'isnl. 

2nd 


(C-7) 


For  cases  other  than  uniform  noise  and  element  power  patterns,  the  reader 
is  referred  to  Table  C-l.  The  required  integration  is  often  difficult  to  attain,  there¬ 
fore  the  resulting  expression  is  written  out  for  each  case  in  Table  C-l. 

Corresponding  to  the  linear  array  case,  a  similar  expression  may  be  defined 
as  the  generalized  correlation  coefficient  for  a  planar  array. 


it 


n,m,  i,  j 


where 


K 


J  °°  j2*  f  1(9, v)  G (9,0)  W(f)  ej  ^  sin  9  d  9  d  0  d  f 
-°°  0  0 


(C-8) 


2n 

\ 


sin  S  [d  (n-i)  cos  0  +  D  (m-j)  sin  0  ] 


(C-9) 


and  K  is  a  normalizing  factor  such  that  Rnj  mj  ^  j  =  J.  for  i  =  n  and  j  =  m.  1(0, 0), 

G(p,  9)  ana  W(f)  are  defined  as  in  equation  (C-l).  The  main  difference  between  equa¬ 
tion  (C-8)  and  equation  (C-l)  is  that  for  the  two  dimensional  case,  the  array  geometry 
must  be  considered. 


Corresponding  to  the  derivations  of  specific  correlation  coefficients  for  the 
linear  array  case,  it  is  possible  to  illustrate  a  mathematical  derivation  of  two  speci¬ 
fic  correlation  coefficients.  These  are  (C-l),  a  planar  array  with  isotropic  radiators, 
a  uniform  spatial  noise  source,  and  narrow  bandwidth;  and  (C-2),  the  same  conditions 
but  with  a  flat  spectral  noise  bandwidth. 


For  the  first  case  we  have: 


G(0,p)  =  1 
I  (0,0)  =  1 
W(f)  =  6  (f  -  fQ) 


(C-10) 


111 


The  correlation  coefficient  is 


R  .  .  =  K  /  /  e 


n,  m,  1,  j 


.  2rr 

2*  n  J  ■>“  sin  e  (d  cos  0  +  d  sin0) 

f  e  *  X  y  sin  0  d  6  d  0 


0  0 


(C-ll) 


where 


d  =  d  -  d. 
x  n  i 


d  =  D  -  D. 

y  m  j 


By  introducing  a  change  of  variable  in  (C-ll)  and  simplifying,  we  obtain: 

_  »,  rn  .  a  ,  .  r2n  j  Z  sin  9  cos  <0-a)  ,  t 

R  .  .  =  K  J  sin  9  d  9  J  e  J  '  '  d0 

n.m.i.j  n  « 


(C-12) 


Fur tlier  simplification  is  possible  if  we  use  the  identity  for  the  nth  order  Bessel 
function  of  the  first  kind  given  by 

2tt 

r  _  j”  t  j  Z  sin  9  cos  0  jn0  . 

J  (Z  sin  9)  s  -J~-  J  e  eJ  v  d0 

n  0 


(C-13) 


Equation  (C-12)  becomes 


R  .  .  =  K  /  2?r  J  (Z  sin0)  sin  0  d  9 


n.m.i.j  '  o 


(C-14) 


or  in  terms  of  original  variables, 


R  .  .  =  K  fn  2rr  J  (¥  sin  9  V  d  2  +  d  2  )sin  0  d  0 
n,  m,  i,  j  o  \  A  x  y  / 


(C-15) 


2tt  /  2  2 

If  Z  =  -r-  v  d  +d  is  small,  we  write  the  Maclaurin  expansion  of 
a  "  x  y 


J  as 
o 


J  (Z  sin 9)  =  ,.,.2  -  ... .2  /0.2 


2  4 

Z  sin  0)  (Z  sin  0) 

...2  ,„,2  +  .2  ~T 


(Z  sin  0) 


(II)  (2)  (21)  (2)  (31  )2  (2)6 


+  .  .  . 


(C-16) 

Evaluating  the  integral  in  (C-15)  by  integrating,  term  by  term,  and  imposing 
on  0  the  limits  0  and  rr,  we  obtain 
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Sin  (r0  Vdx2  +  dy2) 


(C-17) 


R 

n, 


m, 


i.J 


where  the  normalizing  factor,  K,  has  been  applied  so  that  R  .  .  =  1  for  d  = 

n,  m,  i,  j  x 

d^  =  O.  This  equation  is  analogous  to  the  one  for  the  linear  case  (equation  C-7). 

Equation  (C-17)  can  be  extended  to  include  a  fiat  spectral  noise  bandwidth. 
The  correlation  coefficient  becomes 


R 


n,m,  i,  j 


00 

J  W(f) 
-00 


Sin  (ili-  Yd2  ,  J2) 

\  C  '  X  V  / 


2  7rf 
c 


,  2  .  ,  2 

i  +  d 

x  y 


(C-18) 


or 


R 


n,  m,  i,  j 


-p'rH2 

o 


Si 


o  ' 


2n  "  1 '  2  +  dy2  (1  +  P/2) 


(C-19) 


Si 


r  Vdx2  +  d/ 


where  P  is  the  bandwidth  ratio  and  \q  is  the  center  frequency  wavelength. 

Table  C-l  and  Figures  C-l  through  C  -16  give  the  mathematical  form  and 
graphical  representation  of  the  correlation  coefficients  corresponding  to  various 
conditions  of  antenna  element  function,  noise  distribution,  and  signal  bandwith. 
These  examples  are  not  intended  to  be  all  inclusive,  but  representative  of  the  type 
of  variations  which  may  be  expected  in  a  wide  range  of  array  applications. 
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TABLE  C-l.  ARRAY  CONDITIONS  AND  CORRELATION  FUNCTION  (Continued) 


f 


TABLE  C-l.  ARRAY  CONDITIONS  AND  CORRELATION  FUNCTION  (Continued) 
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TABLE  C-l.  ARRAY  CONDITIONS  AND  CC-il'FIATION  FUNCTION  (Continued) 
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TABLE  C-l.  ARRAY  CONDITIONS  AND  CORRELATION  FUNCTION  (Continued) 


TABLE  C-l.  ARRAY  CONDITIONS  A  ND  CORRELATION  FUNCTION  (Continued) 
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i 
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Figure  C-: 
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Figure  C-3.  Correlation  Coefficient 
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Figure  C-5.  Correlation  Coefficient  R  (d) 
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Figure  C~7.  Correlation  Coefficient  R  (d) 
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rrelation  Coefficient  R  (d) 
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Figure  C-9.  Correlation  Coefficient  R  (d) 
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Figure  C-10.  Correlation  Coefficient  R  (d) 
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Correlation  Coefficient  R 
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Figure  C-12 


Figure  C-13 
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igure  C-14 


Figure  C-15.  Correlation  Coefficient  R  (d) 
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Figure  C-16.  Correlation  Coefficient  R  (d) 
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This  study  was  initiated  to  determine  the  effects  of  antenna  element  weighting  (ampli¬ 
tude  control)  and  spacing  on  the  signal-to-noise  ratio  of  wideband  systems.  The 
initial  analysis  Is  directed  toward  the  linear  array  operating  in  an  environment  of 
white  noise,  uniformly  distributed  in  space.  These  confinements  are  later  removed  to 
illustrate  application  to  planar  arrays,  and  non-uniform  noise  conditions. 

It  is  shown  that  the  signal-to-noise  ratio  of  an  antenna  system  operating  in  a  uniform 
noise  environment  is  optimized  when  the  array  directivity  is  maximum.  Optimization 
involves  the  application  of  a  specific  illumination  function  which  can  be  achieved  by 
element  weighting  element  spacing  or  a  combination. 

The  application  of  wideband  signals  modifies  the  array  pattern  in  the  time  space  domain. 
Although  the  nature  of  the  modification  is  dependent  on  waveform,  the  sidelobe  displace¬ 
ment  in  time  and  space  is  a  function  of  the  number  of  elements  comprising  the  array. 

An  iterative  technique  of  matnematical  analysis  is  developed  which  permits  solutions  to 
large  array  problems  where  matrix  inversion  is  difficult.  Several  examples  of  possible 
applications  of  the  analysis  are  presented  including  a  spectral  detection  technique,  an 
approach  to  secure  communications,  and  a  concept  of  array  adaptation  to  a  non-uniform 
noise  environment. 
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